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In this article we review the basic formulation of light-front held theory and light- 
front phenomena in strong interaction. We also explore various approaches to the 
^^O ■ understanding of these phenomena and the associated problems of hadronic bound 

states based on QCD (quantum chromodynamics) on the light-front. 

' PACS. ll.lO.Ef - Lagrangian and Hamiltonian approach. 

t;^- ! PACS. ll.lO.Gh - Renormahzation. 

(N ■ PACS. 11.10. St - Bound state problems. 



I. INTRODUCTION: A HISTORICAL DEVELOPMENT OF 
LIGHT-FRONT DYNAMICS AND SOME COMMENTS 



D , In this article, I will try to address some aspects of nonperturbative QCD. The theory 



of nonperturbative QCD (at the hadronic scale) has indeed not been well defined at the 
present time. Firstly, how much we can understand about the real theory of the strong 
^ \ interaction from the canonical QCD structure is totally unknown. Secondly, there are in 

a I the literature many models or effective theories that are inspired by QCD, and yet they 

cannot be derived directly from QCD, although they are more or less successful in describing 
hadronic phenomenology. Undoubtedly, to discuss the problems of nonperturbative QCD is 
not a simple task. Yet the fact that little progress has been made for the past twenty years 
on the study of nonperturbative QCD (i.e., the explicit solution of the dynamics of quarks 
and gluons in the strong coupling regions) may force us to think about whether one is able 
to extend the traditional approach of the perturbative field theory to the description of a 
strong coupling theory. The covariant perturbative framework lacks the simple Schrodinger 

^ Based on the lectures presented in "the Second Workshop on Particle Physics Phenomenology", May 

19-21, 1994, Renting, Taiwan. 
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picture of quantum mechanics, which may be the obstacle to further development of a 
nonperturbative covariant framework. In this article I am going to introduce a "new" 
dynamical framework, namely light-front dynamics, which can overcome the above obstacle 
and has many advantages for exploring nonperturbative field theory^ that may point to a 
new direction in the investigation of the relativistic dynamics of the strong interaction and 
nonperturbative QCD in the future. 

In field theory textbooks light-front dynamics may not be a familiar term for most of 
us. However, this topic is actually not new at all. In every step of the development of strong 
interaction theory light-front dynamics has indeed played a crucial, if not the essential, role 
for our understanding of hadronic physics. Therefore, before I list the main problems that 
I will discuss in this paper, I would like to give a brief review of the historical development 
of light-front dynamics and the role it played in the study of the strong interaction and 
QCD over the past forty-five years since its invention. However, I must emphasize that 
all I discuss in this article is only a part of the problems in the strong interaction and in 
light-front dynamics that currently interests me. A complete description of the theory of 
the strong interaction and the light-front field theory is certainly beyond my knowledge. 
If the reader is interested in this topic, he or she may find a relatively complete list of 
references in [1]. 

Light-front dynamics, simply speaking, is a description of the evolution of a relativistic 
system along a light-front direction. 

The light-front time-space variables, defined as = t it x where x is one of the 
three components of the space variables, were indeed invented a long time ago, even before 
the birth of modern physics. They have been widely used in solving various problems 
in mathematical physics for the last hundred years. A typical example that one might 
remember from an undergraduate course is how to solve the wave equation (for simplicity, 
we only consider one-time and one-space variable for illustration here): 

where v is the velocity which is a constant. It is well-known that, to solve the above 
equation, it is most convenient to introduce the new variables: ^ = t + x/v and ry = t — x/v. 
With these variables, the wave equation is reduced to ( ggg^ ) ^ = and its solution must 
have the form ^'(t, x) = f{t -\-x/v) + g{t — x/v). Physically, this solution indicates that the 
system moves along the line t = zizx/v with the velocity v in the t — x space. The light-front 
variables x^ correspond to the case where the velocity v = c (the speed of light) with the 
unis c = 1. In such a case, we say that the system moves on the null-plane. 

2 One may note that this framework also has many disadvantages compared to the traditional covariant 
framework if we apply it to perturbation theory, as we will see later. 
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The physical imphcation of the use of the hght-front coordinates can indeed be mani- 
fested from the solution of the above wave equation. First, we consider a small v (compared 
to the speed of light) in Eq. (1.1) which, physically, corresponds to nonrelativistic dynam- 
ics. In this situation, the system moves along the line x = ^vt which does not deviate from 
the line x = very much (see Fig. la). As a result, it may be convenient to describe the 
evolution of the system in terms of the usual time variable t as the time direction. Such 
a choice of the time direction is the one people are most familiar with because this is the 
world wc live in now. However, iiv c (corresponding to relativistic dynamics), the system 
indeed moves near the line x^ = t + x = or x~ = t — x = hi the t — x plane (see Fig. 
lb). In this case, if we still use the variable t as the time direction along which the system 
evolves, the phenomena we measure in laboratories must behave very differently from the 
picture of nonrelativistic dynamics that we intuitively think about. Yet, if we choose a new 
"time" direction, namely in terms of the light- front variable x~^ (or x~), structurally the 
time evolution of a relativistic system may look very similar to a nonrelativistic system in 
the ordinary time-space coordinate system. Consequently, with the light-front variables, 
relativistic dynamics may become simple and transparent. 

Forty-live years ago Dirac first introduced the concept of the light-front form for 
relativistic dynamics, which he called the front form. He defined it as follows: the three- 
dimensional surface in space-time formed by a plane wave front advancing with the velocity 
of light, for example x"*" = t+x^ = 0, is called a front, where x^ is the third component in the 
three-dimensional space. For a dynamical theory, a form of dynamics which is associated 
with the sub-group of the inhomogeneous Lorentz group that leaves the front invariant may 
appropriately be called the front form [2]. The front form is now commonly called the 
light-front framework or simply the light-front. In the literature, it is sometimes also called 
the null-plane. 

Dirac's attempt to introduce the light-front form was to look for the most convenient 
Hamiltonian formulation for a relativistic system. A relativistic system is required to satisfy 
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FIG. 1. The illustration of light-front dynamics. 

two principles: it is invariant under infinitesimal inhomogeneous Lorentz transformations 
and it is a Hamiltonian system. The first requirement arises from the relativistic principle 
and the second is given by quantum mechanics. Then a dynamical system is characterized 
by the ten generators, P'^, M^" , of the Poincare group, which are determined by the Hamil- 
ton equations (Heisenberg equations in quantum mechanics). Based on this philosophy, 
Dirac found that there are three possible forms for relativistic dynamics: i). The instant 
form of dynamics: the dynamical variables referring to physical conditions at some instant 
of time.^ In this case the subgroup that leaves an instant invariant are the translations and 
rotations while the dynamics is determined by and the boost operators M*^*. According 
to Dirac, the instant form has the advantage of being the one people are most familiar with, 
and, that is all. ii). The point form of dynamics: a form of dynamics associated with the 
subgroup of the inhomogeneous Lorentz transformations that leaves a point, say the origin 
= 0, invariant. In this case, the homogeneous Lorentz transformation is simple (kine- 
matic) but all the four-vector momenta P^^ are dynamically dependent. The advantage of 
the point form is perhaps the dynamical separation between P^^ and M^^ . iii). The front 
form of dynamics that has been defined in the previous paragraph. The front form has the 
advantage that of the ten generators of the Lorentz group, only three (instead of four in the 
other forms) are dynamically dependent (as we will discuss later in detail). Furthermore, 
there is no square root for the Hamiltonian, which might simplify the dynamical structure 
since a square root relation between the energy and the momentum cannot provide a simple 
picture of the bound state Schrodinger equation. 

At the time Dirac proposed the light-front Hamiltonian dynamics modern quantum 

^An instant is a flat three-dimensional surface containing only directions which lie outside the light-front 
for example, i = 0. 
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field theory underwent a revolutionary change, namely a covariant perturbative field theory 
was established [3-5]. This formally beautiful and practically useful theory has been used 
so successfully in the description of all the known physical phenomena in quantum electro- 
dynamics (QED). At that time, Dirac's ideas about the front form of relativistic dynamics 
were completely ignored. This should not surprise us because, if we can formulate and 
compute everything in a covariant form, we never need to dynamically solve the Poincare 
algebra. This was perhaps the main reason why the above mentioned ideas of Dirac have 
not recieved enough attention since their invention. Feynman's covariant formulation has 
dominated the investigation of field theory for the past forty-five years. Undoubtedly it has 
so many advantages for our understanding of the microscopic processes with such a novel 
picture and provides a tremendous simplification of the practical computations. However, 
one might recall, that, since it was invented covariant field theory has achieved little progress 
for nonperturbative theory, specially for the relativistic bound state problem. More pre- 
cisely speaking, the quantum field theory we use today may be better to be explicitly called 
the perturbative field theory. Nonperturbative theory requires a complete description of 
the whole physical state space. Yet, in the covariant form, in order to keep the covariance, 
one has to work in a state space that is much bigger than the original physical state space. 
This makes the nonperturbative picture physically unclear and computationally much more 
complicated within the framework of the covariant form. 

Note that, even for perturbation theory, the covariant formalism involves too many 
unphysical processes in practical computations in order to keep its covariant form. A 
typical example from the early days of field theory is the Z-graph (see Fig. 2) which occurs 
in the amplitude of various fermion scattering processes. This diagram contains a state that 
violates the Pauli principle that cannot be explicitly excluded in a covariant calculation. 
In principle, the contribution of this unphysical state is cancelled by the contribution of 
the same unphysical state involved in the vacuum diagram which is, however, factorized 
out in the calculation. In order to avoid such subtle problems, Weinberg developed a new 
formulation for perturbative field theory in the infinite momentum frame in 1966 [7], which 
is equivalent to the field theory quantized on the light-front.^ In the infinite momentum (as 
well as in the light-front) perturbation theory, all diagrams that involve particle production 
from the vacuum must vanish, thereby the Z-graphs are automatically excluded.^ 

In fact, the development of perturbative field theory in the infinite momentum frame 
given by Weinberg was motivated very much by the works of Fubini and Furlan [8] and 
by Dashen and Gell-Mann [9] on the so-called fixed mass sum rules derived from current 

*It must be pointed out that the infinite momentum frame itself is not the same as the Ught-front. The 
former is a special Lorentz frame while the later is a coordinate system of space-time in which all Lorentz 
frames, from the rest frame to the infinite momentum frame, can be described. 

^Strictly speaking, this statement is only conditionally true which we will discuss later. 
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algebras. From the current algebras of the vector and axial vector currents, one can easily 
derive many sum rules for various current matrix elements that can be measured in ex- 
periments. However, in an arbitrary frame, these sum rules may mix the space-like and 
time-like processes and depend on the intermediate state energies so that the momentum 
transfer in these matrix elements in the sum rules varies term by term. As a result, these 
sum rules are difficult to be used to test theory with experiments. Fubini and Furlan found 
that if one boosts the target into the infinite momentum frame, the above difficulty is solved 
and the sum rules are reduced to those with a fixed momentum transfer (g^), namely the 
fixed-mass sum rules. 



Meanwhile, the sum rules become very useful if they are approximately satisfied 
by some set of single particle and resonant intermediate states that provide us with a 
large number of algebraic relations among coupling constants and form factors. In the 
infinite momentum frame the nonzero matrix elements of currents should span an SU{6) 
representation, as was noticed by Dashcn and Gell-Mann [9]. Thus, the main contribution 
to the sum rules comes mostly from the single particles and resonance hadronic states that 
are classified by the SU{6) spectrum generating algebra (the current quark model). 

After the above early application to the study of the strong interaction, the most 
exciting progress made for light-front dynamics was perhaps the parton phenomena ob- 
served in deep inelastic scattering in the very late 60's. The parton picture, as everyone 
knows, was introduced by Fcynman [10,11] for interpreting the scaling behavior discovered 
by Bjorken [12] from the deep inelastic scattering experiments. Although most of us will 
immediately think of the partonic picture in the infinite momentum frame, when we talk 
about parton phenomena, it is in fact much easier to see its theoretical underpinning on the 
light-front. The cross section for deep inelastic scattering is proportional to the so-called 
hadronic tensor which is a Fourier transformation of the hadronic matrix element of the 
commutator of two currents: 



In the rest (Lab) frame with high energy and high momentum transfer where the transfer 



FIG. 2. The Z-diagram. 




(1.2) 
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momentum is assumed to be along the z-direction, it is easy to check that the dominant 
contribution in the above integral comes from the region = t + z ^ 0. In other words, 
in high momentum processes, the main contribution to the hadronic tensor comes from the 
light-front commutator of two currents. On the other hand, causality, = x'^x~ — x^>0, 
shows that the light-front commutator is not vanishing only if x± = 0. Obviously the 
commutator of the two currents on the light-front can only depend on x~ , the light-front 
longitudinal variable, from which one can easily derive the Bjorkcn scaling that the structure 
functions extracted from the hadronic tensor are independent of the hadronic mass scale. 
Gross and Treiman further proved that the dominant contribution to the hadronic tensor 
is from the most singular part of the light-front commutators which is mainly given by 
the free theory [13]. This provides a theoretical foundation for Feynman's parton model. 
Therefore one usually refers to the parton phenomena the light-cone dynamic of the strong 
interaction. 

The manifestation of light-front dynamics from the parton picture in deep inelastic 
scattering led to many attempts to reformulate field theory on the light-front in the very 
early 70's. Typical works are by Drell, Levy and Yan [14-16] and by Bjorken, Kogut and 
Soper [17], who attempt to derive the parton structure directly from field theory [18]. 
Meanwhile, Chang and Ma began to study the canonical structure of field theory on the 
light-front [19]; Leutwyler et al. attempted to develop a light- front quantization approach 
to field theory [20]; Kogut and Soper reformulated QED on the light- front [21]; Rohrlich 
et al. explored the relation of the initial value problem with boundary conditions on the 
light-front for QED [22,23], and Chang and Yan extensively discussed various field theories 
quantized on the light-front using Schwinger's action principle and compared the light-front 
noncovariant perturbative structure to the traditional covariant formulation [24-27]. Also 
there were many other concerns in the light-front field theory at that time (for examples, 
see [28,29]). These investigations provide a basis for conventional canonical light-front field 
theory. However, the nonperturbative aspects of field theory on the light-front had not 
been addressed at that time. 

On the other hand, the parton phenomena also led to further extensive exploration 
of current algebras on the light-front and the associated sum rules. For example, the fixed- 
mass sum rules in the infinite momentum frame can be easily and naturally obtained from 
the light-front current algebras, which were provided by Jackiw and his collaborators [30,31]. 
It was also found that with the light-front current algebras some anomaly structures, such 
as Schwinger terms, may be easy to address [32,27]. The most important development for 
strong interaction along this line is the development of a fundamental theory for the strongly 
interacting basic constituents (quarks and gluons), i.e., QCD. Now everybody believes that 
QCD is the only acceptable fundamental theory for the strong interaction, but maybe not 
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everyone knows that QCD was initially formulated by Gell-Mann and Fretsch [35] on the 
light-front. At the time many particle theorists focused on the field theory interpretation 
of the parton phenomena, Gell-Mann and many others still concentrated on looking for a 
possible fundamental theory for the quark model which is very successful in the description 
of the classification of hadronic states. In the current quark model, an essential equation 
to determine the dynamics of the model is the quark equation of motion on the light-front. 
Gell-Mann and Frctszh wrote down the quark equation of motion exactly, based on SU (3) 
gauge theory on the light-front with the light-front gauge, in their paper of 1972. 

Of course, QCD is accepted as the fundamental theory of strong interactions due 
mainly to the discovery of its asymptotic freedom behavior by Gross and Wilczek and by 
Politzer [36,37], which provides the theoretical basis for the parton phenomena in high 
energy physics. Asymptotic freedom allows us to use perturbation theory to do a QCD 
calculation for physical processes near the light cone ~ (the common statement that 
perturbative QCD is for short distance dynamics is indeed a misleading comment). Ob- 
viously, the natural framework for exploring the parton phenomena with QCD is the for- 
mulation of QCD on the light-front with the light-front gauge. The remarkable properties 
of the light-front formulation are that the conventional vector and axial vector currents of 
quarks remain gauge invariant in QCD, and the concept of parton distribution functions 
naturally occurs [38]. 

The development of a light-front formulation for QCD is as follows. Tomboulis 
first quantized Yang-Mills theory on the light-front in 1973 [39]. Cornwell and Crewther 
discussed possible singularites arising from the choice of light-front gauge [40,41]. The 
first QCD calculation with the light-front gauge was given by Gross and Wilczek in 1974 
[42], who used it to check the gauge dependence of the one-loop anomalous dimension 
and the P function. The canonical quantization of QCD on the light-front was provided 
by Casher [43] and by Bardeen et al. [44,45] in 1976. Thorn addressed the problems of 
asymptotic freedom and the confinement mechanism within light-front QCD in the 
late 70's [46]. Lepage and Brodsky developed the Hamiltonian formulation of light-front 
QCD and explored the asymptotic behavior of hadrons and the perturbation expansion of 
various exclusive processes [47-49] . Franke et al. first addressed the problem of light-front 
gauge fixing and the associated nontrivial QCD structures [50]. Glazek discussed the quark 
condensate in light-front QCD [51]. Meanwhile, due to the scaling evolution of the hadronic 
structure functions suggested by Altarelli and Parisi [52] , a systematic development of the 
factorization theorem for various inclusive processes and a practical calculation scheme of 
the various anomalous dimensions within light-front QCD were also provided in the early 
80's [53-55]. Note that the factorization theorem for deep inelastic scattering was initially 
proven by the use of operator product expansions which are only valid near the light-cone 
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[56]. In the mid 80's, there were also a lot of discussions about the trouble caused by the 
light-front gauge singularity and its possible solutions in perturbation theory [57-60]. In 
conclusion, due to the light-cone behavior of high energy dynamics, perturbative light-front 
QCD, as the foundation of the parton picture, has been widely used in various investigations 
of high energy processes. 

Here I shall make some comments. Because of the two conventional requirements of 
field theory, namely, covariance and unitarity, the development of non-abelian gauge theory 
with gauge fixing is not at all simple. In the usual covariant path integral quantization, 
the unpliysical gauge components violate unitarity and therefore one has to introduce ghost 
fields to cancel the unphysical states in the theory [61]. A covariant formalism has several 
advantages for many aspects, yet the physical picture becomes obscure once unphysical 
particles are introduced. On the other hand, in canonical quantization one has to choose a 
physical gauge, such as the axial gauge or light-front gauge, where unitarity is automatically 
satisfied but covariance is no longer manifest. Although we know that the natural picture 
of QCD for parton phenomena is based on the light-front gauge, many still prefer to use 
the covariant formulation. The main complaint against the light-front gauge is that the 
usual prescription (principal value prescription) for the gauge singularity prohibits any 
continuation to Euclidean space (Wick rotation) and hence power counting for Feynman 
loop integrals. As a result, non-local counterterms have to be introduced, which break the 
multiplicative renormalizability [41,59]. 

However, even in covariant perturbation theory, formal multiplicative renormalizabil- 
ity may not be significantly useful in the calculation of hard scattering coefficients. This 
is because gauge invariant composite operators in the operator product expansion are not 
multiplicatively renormalizable [62], and the gauge variant operators will induce gauge in- 
variant counterterms that do contribute to gauge invariant matrix elements. Although it 
was believed that the mixing problem had been solved [63], the long-standing discrepancy 
between the covariant and light-front gauge calculations for the second-order anomalous 
dimension of the gluon composite operators has recently renewed attention to the renor- 
malizability of QCD [64,65]. In fact, the mixing of gauge invariant operators is an intrinsic 
property that exists in any covariant formulation of QCD. The above mixing problem of 
gauge invariant operators apparently disappears in the physical gauges and perturbative 
calculations become straightforward. However, a new spurious mixing of ultraviolet and in- 
frared divergences associated with the gauge singularity does occur in the light-front gauge 
using the principal value prescription. One has to be careful to see whether or not the 
spurious mixing affects physical quantities. In the light-front gauge, it has been shown that 
the most severe infrared divergences caused by the A'^ = gauge choice with the principal 
value prescription are indeed cancelled in the gauge invariant sectors [53]. Thus, in the 
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perturbation region, the light-front gauge singularity in QCD may not be a severe problem. 

Although QCD manifests asymptotic freedom, the main feature of QCD as the fun- 
damental theory of the strong interaction should lie in its nonperturbative behavior. In 
the past twenty years, there has been tremendous work on nonperturbative QCD, and 
some nontrivial structures, such as the instanton [66], were revealed. Unfortunately, the 
real picture of nonperturbative QCD that leads to quark confinement and chiral symmetry 
breaking, which are believed to be the two essential mechanisms for form hadron formation 
from QCD, is still lacking. Roughly speaking, all methods that we use today for explor- 
ing nonperturbative QCD were mainly developed by Wilson in the late 60's and the early 
70's. These are the lattice gauge theory [67], operator product expansions [56] and the 
renormalization group approach [68,69]. The lattice gauge theory is the only computable 
approach in the market for the nonperturbative study of QCD. In the past one and a half 
decades, extensive work on lattice gauge theory has led to significant progress in many 
research avenues; however, accurate information on the bound state of light quarks is still 
not available [70]. On the other hand, a phcnomcnological approach of nonperturbative 
QCD, namely QCD sum rule [71], has become a reliable approach for various hadronic 
structures. The theoretical basis of QCD sum rule is the operator product expansion which 
separates the perturbative and nonperturbative contributions of QCD in hadrons near the 
light-cone. In fact, what one has learned from QCD sum rules about QCD is nothing more 
than the perturbative part of QCD. The nonperturbative aspect in QCD sum rules is still 
purely phenomenological, including the parameterization of the quark and gluon condensate 
(which must be directly calculated from QCD for a full understanding of nonperturbative 
QCD but so far no one knows how to compute them). 

Succinctly speaking, the aim of nonperturbative QCD is to find the hadronic bound 
states directly from QCD. However, the only truly successful approach to bound states 
in field theory has been Quantum Electrodynamics (QED), with its combination of non- 
relativistic quantum mechanics to handle bound states and perturbation theory to handle 
relativistic effects [72-75]. Lattice Gauge Theory is maturing but has yet to rival QED's 
comprehensive success. There arc three barriers which prohibit an approach to Quantum 
Chromodynamics (QCD) that is analogous to QED. They are: (1) the unlimited growth of 
the running coupling constant g in the infrared region, which invalidates perturbation the- 
ory; (2) confinement, which requires potentials that diverge at long distances as opposed to 
the Coulombic potentials of perturbation theory; (3) spontaneous chiral symmetry break- 
ing, which does not occur in perturbation theory. The last two barriers are essentially 
related to the nonperturbative structure of the QCD vacuum. 

In contrast to the gloomy picture of the strong interaction in QCD, however, is that 
of the Constituent Quark Model (CQM), the most successful model in the description of 
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hadronic dynamics. In CQM, only a minimum number of constituents required by the 
symmetries is used to build each hadron [76]. Moreover, Zwcig's rule leaves little role 
for the production of extra constituents [77]. Instead, the rearrangement of pre-existing 
constituents dominates the physics [78-80]. Yet the CQM has never been reconciled with 
QCD — not even qualitatively. Note though that QCD was initially proposed on the light- 
front in order to develop a fundamental theory for the quark model [35]. In 1989, Wilson 
suggested that QCD on the light-front may provide a natural framework to solve QCD 
nonpcrturbativcly for hadrons because the use of the light-front coordinates may allow us 
to develop the same techniques that solve QED for solving the physics of CQM [81]. To 
do so, he proposed a new idea for constructing the low-energy QCD theory for hadrons 
based on light-front power counting which differs from what we are familiar with in the 
equal-time framework. Combined with the two-component formulation of the light-front 
QCD [82] and a new scheme for Hamiltonian renormalization [83,84], it seems to provide a 
newer avenue for exploring nonperturbative QCD dynamics and hadrons than that used for 
solving QED. This has been presented in a recent paper by Wilson, Walhout, Harindranath, 
Zhang, Perry, and Glazek [85]. 

At this moment, I cannot predict whether this new approach will give the final 
solution for QCD or whether this approach will be successful for a strong coupling theory. 
Light-front dynamics may shorten the distance from our hopes for QCD to a real picture 
of the hadronic world. What I can comment on is that in the past twenty years various 
phenomenological understandings of nonperturbative QCD, such as the constituent quark 
model for the heavy quark system, the heavy quark symmetry for the heavy-light quark 
system and the effective chiral theory for the light quark system, could guide us in the right 
direction. Optimistically, we should believe that any new attempt inspired by a successful 
understanding of phenomenology always contains a possibility of providing a real answer 
to the unsolved problems. 

After this long historical overview of the development of light-front dynamics and 
light-front QCD in terms of my personal point of view, I now give an outline of this review. 
In the next section, I will present a brief description of the basic structure of the light- 
front form and the free field theory on the light-front. Then in Sec. Ill, I will discuss 
the light-front phenomena from which I hope a picture is provided of what we really need 
to understand from QCD on the light-front for strong interaction phenomena, i.e. the 
light-front bound states for hadrons. In Sec. IV, I will discuss the general structure of 
light-front bound states and their advantages in describing hadronic physics, and some 
phenomenological light-front wave functions for hadrons that have been used for successfully 
describing many hadronic properties. In Sec. V, as a basis for exploring the light-front 
dynamics of QCD, canonical light-front QCD and its various properties and problems are 
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discussed. In Sees. VI and VII, I will introduce new ideas that were originally proposed by 
Wilson and his collaborators for finding light-front QCD bound states for the hadrons, and 
I will also discuss possible practical schemes to realize these ideas and thereby to provide a 
possible real theory of QCD for hadrons. 



II. LIGHT-FRONT FORM 



II-l. Lorentz Transformations on the Light-Front 

One of the useful properties in light-front formulation of relativistic dynamics is its 
convenient Lorentz transformation structure [2]. In this section, we shall introduce some 
basic properties of the Lorentz transformation on the light-front. 

The light-front coordinates of space-time — , — , — ) are defined as 

x±=x°±x^ x\ = x\ (z = l,2), (2.1) 

where we choose x"*" as the "time" direction along which the system evolves. Thus and 
x_L naturally become the longitudinal and transverse coordinates. All other 4-vectors in 
space-time are defined in the same form as Eq. (2.1). The inner product of any two four- 
vectors is given by a^b^^ = ^{a^h~ + a~h^) — a_\_ ■ b±, where the metric tensor for Eq. (2.1) 
is = g = 1/2, g^^ = —6^^ and the other components are zero. The light-front time 
and space derivatives are d~ = 2-^|qr, = 2-^^, 5* = = — -^r (see the definition of 
the metric tensor), and the four-dimensional volume element d'^x = ^dx~^dx~(Px±. 

1. Poincare algebra. A relativistic dynamical system must be Lorentz invariant. The 
inhomogeneous Lorentz transformation is generated by the Poincare algebra: 
[p/^, pf] = 0, [M^"',PP] = ii-gt^pP" + g^ppf) , 

(2.2) 

[MP", MP''] = ii-gPPM"" + g^PMP" - gP" MP" + g'^'MPP) , 

where PP is the four-vector momentum, and MP" describes the rotational and boost 
transformations. In the instant form, the rotational and boost operators are given as 
M*^' = eijkJ'' and = K\ respectively, with 

[J^ = ieijkJ'' , [J\K^] = ie^jkK'' , [K\ K^] = -ieijkK'' . (2.3) 

In field theory these generators are constructed from the energy momentum density 
QP" . For a given Lagrangian C{(j)i), the energy-momentum tensor in the instant form is 
defined byf] 

d£ 

= ^r^d-'^^-gP'-C. (2.4) 

®If C contains gauge fields, we must add to O'"' an additional term dp{F'^'^ A") to ensure that it is gauge 
invariant, symmetric and traceless. 
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Then, 

P^ = j d^xQ'^" , M'^'^ = j (fxlx"®'^" - x'^e^^l . (2.5) 

In light-front coordinates, the four-vector momentum operators and the boost and rota- 
tional operators are defined by: 

p± = ±P^ , Pi = P' {i = 1,2), 

= K^ + J^, E^ = K^-J\ K^, (2.6) 

In terms of the energy-momentum density, = —^M'^ , E'' = and = M^^, 

F"- = M-\ where 

P^ = \J dx-d^x±e+^' , ^4^" = ^ J dx-d^x±[x^'e+'' - x'^e+'^j . (2.7) 

Physically, the light-front transverse boost (-E") is a combination of a pure Lorentz boost 
and a usual rotation, and so is the light- front transverse rotation (F^). 

2. Boost transformation. The most useful property of the Lorentz transformation on 
the light-front is the kinematics of its boost transformation. The subgroup that leaves the 
light-front x"*" = invariant is the translations (P"*", P^, P^) and the Lorentz transformations 
generated by J^, K^, E^. A boost transformation, 

L{p) = eM-iPiE'} exp{-i/33K^} , (2.8) 

transforms the light-front momenta operator as follows: 

L(P)P+L-\p) = P+e* , L(P)P±L-\P) = P±+ P+/3±e^3 _ (2.9) 

Consider a state in the rest frame: |P), where P = (M, 0,0,0) in the usual instant coordi- 
nates. The light-front momentum P"*" = M, Pj_ = 0. Let the transformation parameters 
be f3± = p±/p'^, (3s = ln(p"'"/M), then we have 

L{p)\P) = \P') , P'+ = p+. Pi = , (2.10) 

which boosts the state in the rest frame to the frame with momentum P' = p. The boost 
for a state with an arbitrary momentum A is given by 

L{P)\A) = \A') , A'+ = {P+/M)A+ ; A'^ = A± + {A+/M)P± . (2.11) 

The remaining three generators of the Poincare algebra do not keep the light-front 
plane invariant. The operator P~ transforms the plane x"*" = to the plane x'^ = t while 
P* rotates the plane x'^ = around the surface of the light-cone x"^ = 0. In other words, 
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once we specify the initial data on the hght-front = 0, the three operators, P~ and 
F*, determine the evolution of the system away from x"*" = 0. Obviously, x"*" = r is the 
evolution parameter, namely the light-front time variable. 

II-2. Basic Formulation of Light-Front Field Theory 

In this section we discuss some basic light-front structures in the field theory by 
looking at the canonical light-front quantization for a free scalar, a free fermion and a free 
gauge theory. The motivation is to show how light-front field theory differs from instant 
field theory. 

1. Scalar field. Consider a free scalar field theory, 

-^ = ^(9^^-"^^), (2.12) 

the equation of motion on the light-front is given by 

{d~d+ - dl+m^)^ = Q. (2.13) 

The light-front canonical theory can be obtained by defining the conjugate momentum of 
(j) referred to the light-front time x'^: 

It is worth noting that vr^ is not a light-front time derivative of cj), which implies that the 
light-front conjugate momentum is not an independent dynamical degree of freedom. The 
equal light-front time canonical commutation relation is then given byQ 

['t>{x),T^^{y%+=y+ = (^(a;" - y')5'^{x± - y±) . (2.15) 
If we define 

—jf{x-) = -J^^dx'-e{x--x'-)f{x'-), (2.16) 
where e(x) = 1, 0, —1 for x > 0, = 0, < 0, then 

[cPix), cPiy)U=y+ = -j<x- - y-)6\x^ - y^) . (2.17) 

^Since the conjugate momentum is not an independent dynamical degree of freedom, the hght-front dy- 
namical system is always a constrained Hamiltonian system for which the naive canonical quantization is 
generally not valid, and one should use the Dirac quantization procedure [86] . However, here the con-straints 
can simply be solved so that the naive canonical approach, or more generally the phase space quantization 
approach is good enough to provide the basic commutation relations without any loss of the generality 
[87,88]. 
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The light-front Hamiltonian of the free scalar field is given by 
Hlf = P- = j dx-(fx±{ir^d-(j) - C) 



(2.18) 



A solution of the free scalar field on the light-front in terms of momentum space is 

2(2vr)3[fet] ["(^)^'''' + ' ^ ' (2-19) 

with the commutation relation: 

[a{k),a^{k')] = 2{2Trf k+5{k+ - k'+)5'^{k± - k'^) . (2.20) 

In momentum space the light-front Hamiltonian becomes 

f dk^(Pk\ k'^\ + i ^ ^ , , , , 

where an infinite zero-point energy has been ignored. Thus, the light-front energy k~ for 
the scalar physical state \k) = a'''(A;)|0) is 

HLF\k) = k-\k) ^k- = . (2.22) 

This is the light- front energy dispersion relation, a result of k'^ = k^k^ ~ ^± — 

The motivation for presenting such a detailed discussion of the free scalar theory here 
is to show the difference between the light-front field theory and the instant field theory. 
There are three major differences we can see from the above discussion, i) On the light-front, 
the canonical conjugate momentum is not a dynamical degree of freedom. In other words, 
the number of the dynamical degrees of freedom is reduced by half compared to the instant 
formulation. As a result, the scalar field variables themselves do not commute with each 
other on the equal light-front time plane (see Eq. (2.17)). ii) The energy dispersion relation 
does not involve a square root structure which may help us to formulate the relativistic 
bound states in terms of a simple nonrelativistic-like Schrodinger picture, iii) There is 
a new singularity in the light-front formulation which exists even in the free theory and 
which is hidden in the operator In momentum space it corresponds to the singularity 
at /c"'" = in the form . As we have mentioned in the introduction and as we will discuss 
later in detail the k^ = singularity is the most difficult problem to handle in the light- 
front field theories but it also plays an essential role in determining the nontrivial structure 
of the theory on the light-front. The definition of Eq. (2.16) is to remove this singularity. 
The notation denotes the Fourier transformation of Eq. (2.16): 

1 - 1[ ' I M (2 23) 

[k+] 2[k+ + i€ k+-ie\' ^ ' ^ 



731 



732 



LIGHT-FRONT DYNAMICS AND LIGHT-FRONT QCD 



VOL. 32 



which corresponds to a principle value prescription. Since Eq. (2.16) is an integral it allows 
an arbitrary constant in defining the operator 1/9+ . Different definitions of 1/5"*" can 
lead to different commutation relations of Eq. (2.17) and different prescriptions for pp. 
However, Eq. (2.16) has two advantages: it uniquely determines the initial and boundary 
conditions for the light-front field formulation [22], and it explicitly eliminates the A;+ = 
modes from the theory [88]. The latter property plays an important role in simplifying 
the vacuum structure for light-front bound states. We shall come back to address these 
properties in Sec. V where we shall discuss them in nonabelian gauge theory. Keeping 
these differences in mind, we now discuss the fermion field formulation on the light-front. 

2. Fermion field. For a free massive fermionic field with spin ^ the symmetric 
Lagrangian is 

C = V^(i7- H -m)i} , (2.24) 
which leads to the Dirac equation 

{i-f-d- m)il) = , V'(«7 ■d + m) = Q. (2.25) 
In the canonical procedure the light-front conjugate momentum of and V'^ are defined by 

where From the above expression, we see that if we define 

V'±=A±^, A± = i7V, (2.27) 

where Aj- are the light-front fermionic projectors: A^ = A+, h?_ = A_, A+A_ = 0, then we 
have ijj = il>+ + ijj- and tt^ = yV'+i ^r^t = 1^ other words, the fermion field on the 

light-front can be divided into the up {+) and down (— ) components, while the canonical 
conjugate momentum of ip only relates to the up component. 

The most interesting property for the light-front fermion field is that the Dirac equa- 
tion becomes two coupled equations: 

id-tp+ = {ia± ■ d± + (3m)ip-. , (2.28) 

id+ip- = {ia± ■ d± + f3m)^+ . (2.29) 

The equation of motion for V- turns out to be a constraint equation since it does not 
contain any light-front time derivative. Using the definition of Eq. (2.16), we can express 
V'- in terms of V'-i- and then substitute it into the equation of motion for ip^, obtaining 

id~il;+{x) = — dx'~e{x~ - x'~){-d'i + m'^)'tp+{x+,x'~,x±) . (2.30) 

^4 J —oo 
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or 



{d+d- -81 + m'^)'ilj+ = . 
If we introduce the following 7 matrix representation: 



(2.31) 



7 



7 



_ 



-i 

1 

a' ( 




7 



i 

1 



-ia' 
ia^ 



-a 



(2.32) 



then 



A 



+ 



"10" 


, A_ = 


"00" 







1 



the Dirac spinor is reduced to 



V 



+ 









V 



-^{(T±d± +m)ip 



(2.33) 



(2.34) 



In other words, the fermion field on the light-front can be expressed by a two-component 
spinor ip which satisfies the Klein-Gordon equation, and all 7 matrices can be reduced 
to the Pauli matrix. The relativistic fermion particle can then be described in terms of 
a nonrelativistic spin y particle on the light-front. This is a very useful property of the 
fermion field on the light-front. The canonical commutation relation can be simply obtained, 



5{x - y )5^{x± - y±) . 



(2.35) 



It must be emphasized that the above two-component light-front field still describes 
relativistic spin- 1/2 particles, which are intrinsically different from spin- 1/2 non-relativistic 
particles. In other words, it contains quarks and antiquarks. To see how antiquarks can be 
described in a two-component representation, we need to construct the charge conjugate for 
if. By using the condition Cj'^'^C~^ = it is easy to find that the charge conjugation 
operator in the 7-representation of Eq. (2.32) is 



C = i^l 







la 







Hence, 



,0„/,* 



(2.36) 



(2.37) 



where ry is a phase factor which will be set to unity in this paper. Eq. (2.37) determines the 
definition of fermion and antifermion in a two-component representation on the light-front. 
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The light-front Hamiltonian of the free fermion field is given by 

Hlf = J dx-cFx±{Tr^d-tP - ^) = J dx-cfx±ip^ '^~^]"g"t "^^^ "P- (2-38) 
The solution of the free fermion equation of motion is 

^(x) = X:Xa / ^^^[b{p,X)e-'P' + d^{p,-\)en, P+>0, (2.39) 

where X1/2 = (q)' ^-1/2 = (1) &re the spin eigenstates, and the fcrmionic and antifermionic 
annihilation and creation operators satisfy the following commutation relations: 

[b{p,x),b^ip',x')] = [d{p,x),Sip',x')] 

= 2(27r)3,5A,A'<5(p+ - P+)6\p± - p'J . 
In momentum space, the light-front Hamiltonian is 



(2.40) 



X •' 2(27r)'^ \p+\ ^2.41) 
p+ > 0, 

where again an infinite zero-point energy has been ignored. Thus, the light-front energy p" 
for a fermion physical state \p,X) = b'^{p,X)\Q) is 

2 2 

Hlf\p, X) = p-\p, X)^p- = ^^^^ . (2.42) 

The same result can be obtained for the antifermions. This is the light-front energy disper- 
sion relation for the fermion field which is the same as that for the scalar field. 

3. Gauge field. In the last part of this section, we consider a free gauge field: 

jC = - ^F^.F''^ , F^, = d^A, - d^A^ . (2.43) 
The light-front canonical momentum of A^^ is defined by 

It immediately shows that E'^ = 0, which implies that A~ is not a physically independent 
field variable. With these canonical conjugate momenta, the equations of motion, d^F^^^ = 
0, are reduced to the light-front Maxwell's equations: 

^d+E- -d±-E± = 0, d'E- = d^B^ - d'^B^ , 

(2.45) 

d-E^ = —d+B^ + d'^B- , d-E^ = -—d+B^ - d^B- , 

2 ' 2 ' 
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where E = {E-,E^,E'^), and B = {B' = F'^'^,B^ = F^-,B'^ = F"^). 

Gauge invariance allows us to choose a special gauge. The most convenient gauge 
choice in light-front gauge theory is the light-front gauge 

= A^ + A^ = 0. (2.46) 

With this gauge, the canonical momentum is reduced to E~ = —^d'^A", E^ = —^d'^A^. 
Thus the first light-front Maxwell equation (Gauss Law) turns out to be a constraint equa- 
tion, 

^d+A- = d'A' , (2.47) 

from which the A~ component can be explicitly determined in terms of the physical (trans- 
verse) gauge field. Therefore, a gauge theory on the light-front can be expressed explicitly 
with the physical gauge field A'-: 

Hlf = I dx-cfx^ ^{d'A{f . (2.48) 

The canonical commutation relation for the physical gauge field is the same as that 
of the scalar field: 

[A'{x),A^{yX+=y+ = -^S,je{x- - y-)6'{x^ - y^) . (2.49) 
The solution of the free light-front gauge field in momentum space is 

^'(^) = E/^|^^[^W«'^)e"''' + M, 9+>0, (2.50) 

where e\ is the gauge field polarization vector: e\ = "^(1)0) ^-i = IJT^'^' 
canonical commutation relation for vector boson annihilation and creation operators is: 

[a{q, A), at(g'. A')] = 2i27rfq+5x,x'S{q+ - q'+)S\q^ - q'^) . (2.51) 

In momentum space the light-front Hamiltonian has the same structure as the other fields. 

The light-front energy q" for the physical gauge state \q, A) = b^q, A)|0) is then 

HLF\q,\)=q-\q,\) = (2.53) 

Namely, = 0. 
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It may be worth noting that there is a remarkable property of the hght-front free 
field theory discussed in this section: both fermion and gauge fields are reduced to a two- 
component representation. Later we will see that this is still true even for QCD. In QCD, 
since the gauge fixing = removes the unphysical degrees of freedom, unitarity is 
automatically satisfied and no ghost field is needed, the price to pay being the loss of 
manifest covariance. The choice of the A'^ = gauge promises that QCD involves only 
two-component gauge fields and two-component quark fields in light-front coordinates. In 
early stages of the development of light-front QED, Bjorken, Kogut and Sopcr used the two- 
component formulation to discuss various physical QED processes [17]. In the middle-80's 
a two-component Feynman perturbation theory with the Mandelstam-Leibbrandt prescrip- 
tion [57,58] was developed by Capper et al. [59] and Lee et al. [60] for pure Yang-Mills 
theory. It has been shown that there are many advantages in such a two-component light- 
front field theory. The extension of the two-component theory to the study of light-front 
canonical quantization (i.e. ,t "'"-ordered theory) of QCD has been explored just recently 
[82]. Practically, one will sec that the two-component formulation indeed provides a trans- 
parent physical picture. In order to understand clearly the origin of light-front singularities 
in various physical processes and to simplify practical calculations, all discussions in this 
article are based on the two-component formulation. 

With the above basic formulation of the light-front field theory, we now discuss 
physical phenomena of the strong interaction on the light-front. 

III. LIGHT-FRONT PHENOMENA 

It has been well-known as early as the 60's that the partially conserved axial cur- 
rent (PCAC) provides the dominatant theoretical underpinning for the strong interaction. 
Although QCD has now been commonly accepted as the fundamental strong interaction 
theory, we are still unable to solve hadronic structures from QCD in the low-energy scale 
due to the lack of a practically computable nonperturbative approach for field theory. In 
the past three decades, most hadronic structures at low-energy can only be extracted from 
various phenomenological matrix elements of the currents in the hadronic bound states 
and the sum rules that these matrix elements satisfy by the use of current algebras. In 
other words, after about thirty years of investigation of the strong interaction for quarks, 
it may still be true that of all the theoretical methods which have been invented to de- 
scribe hadrons, only the algebra of conserved and partially conserved quark currents has 
led to precise connections between the underlying fundamental degrees of freedom and the 
observed hadron properties. Typical examples are: i), the low-energy theorem for hadrons 
developed from current algebras in the middle of the 60's which underlies almost all of the 
theoretical understanding of the low-energy strong interaction; ii), the parton dynamics for 
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the deep inelastic scattering invented in the late 60's which has dominated the study of 
high-energy hadronic physics in the last two and a half decades; and iii), the QCD sum 
rules of the quark current correlations for meson and baryon structures investigated in the 
80's, which are currently a main approach of the phenomenological nonperturbative QCD 
descriptions for hadrons. In fact, it is particularly of interest to see that the connections 
between the fundamental degrees of freedom and hadronic observables via currents are nat- 
urally manifested when we use light-front coordinates. Furthermore, the most important 
feature of the above connection in the light-front coordinates is that it directly addresses the 
fundamental interactions between the underlying degrees of freedom (quarks and gluons) in 
terms of physical observables that may further provide an alternative but more transparent 
approach to the study of nonperturbative QCD dynamics for hadrons. We will discuss this 
feature later. In this section I will review some phenomenological descriptions of strong 
interaction on the light-front. 

III-l. Fixed-mass Sum Rules and Light-Front Current Algebra 

The earliest application of light-front dynamics to strong interaction phenomena 
is the derivation of fixed-mass sum rules from current algebra in the infinite momentum 
frame. An introductory discussion can be found in the famous book by Adler and Dashen 
[89] . Without any loss of generality, let us consider the matrix element of the equal-time 
commutator of two vector currents sandwiched between the spin-averaged hadronic states: 



The sum rule can be obtained by inserting a complete set of intermediate states {|n)} into 
the left-hand side of Eq. (3.1): 



(27r)3^(|(/3(p)|j0(0)|n)p-|(/3(p)|j0(0)|n)|2)<53(p + q-p„) = 4/l3(/?), (3.3) 



where the momentum transfer = (p — Pn)^- The difficulty with the above sum rules is 
that depends on the energies of the intermediate states. These matrix elements can be 
measured or estimated from PC AC for certain values of g^, while in the sum over states 
q^ increases without limit as the mass of |n) increases. There is no obvious condition for 
the convergence of this summation. Consequently, its applicability is not clear since, in 
practice, one has to truncate the sum over states. 

The difficulty can actually be removed if we boost the hadronic states to the infinite 
momentum frame, as was first discovered by Fubini and Fur Ian [8]. In the infinite mo- 



(/3(p)|[J+(q), J_(-q)]|/3(p)) = 2(^(p)|J3|;3(p)) , 



(3.1) 




(3.2) 



n 
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mentum frame, |p| oo. Also, for simplification, it is convenient to impose the condition 
p • q = 0. Then 

qO = (v'(p + q)2+^2 _ ^p2 +rn} f I^J^^ Q , (3.4) 

so that q'^ becomes independent of the mass of the intermediate state and is fixed: q'^ = — q^. 
All matrix elements in the sum rule have then a fixed momentum transfer. This leads to 
the so-called fixed-mass sum rule, 

/ duW2{y,q^ = -ci^)=mi3), (3.5) 
where i' = p ■ q, and W2 is defined by 

C^^''ip,q) = i2nfY.if^ip)\JU0)\n){n\ri0mp))5Hp + q-Pn) 

(3.6) 

= -gf^'Wi + p''p''W2 + ql^q^Ws + —{p'^q" + p''q^')Wi 

in which Wj are the structure functions measured in experiments. This is perhaps the first 
application of the infinite momentum frame to the strong interaction phenomena.^ 

However, one may note that in deriving the fixed-mass sum rule, one has in fact used 
the assumption that it is legal to take the limit inside the integral. The validity of the above 
assumption is questionable due to the possible illegitimate interchange of limit and integral. 
It was shown that the fixed mass sum rules are equivalent to the appropriate light-front 
commutators and to the |p| 00 technique [90,91]. One can derive the various fixed mass 
sum rules from the light-front commutators without using the unreliable |p| 00 [31]. 

The light-front current algebra was given by Cornwell and Jackiw [30] based on the 
basic commutation relation for the light-front fermionic field: {^/'+(x), 'i/'+(2/)} = 6{x~ — 
y~)(5^(x_L — y±). For the SU{3) currents, we obtain 

-^5ai,dtd+[e{x- - y-)5\x^ - y^)S] , 

*The above fixed-mass sum rules can also be obtained in the framework of a dispersive formalism with 
an unsubtracted assumption [89]. 



(3.7) 
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= ifabcVc{x)5{x~ - y~)S'^{x± - y±) 

-^{ifabc + dabcKdtHx- - y~)S^{x± - y±)V-{x\y)] 

—diHx- - y-)5\x^ - yA_){V^{x\y) + ie'^ A^Mv)]} ^^'^^ 

- y')5'^{x^ - yi_)'ip{x)-i^-i~ Kabi^iy) - /i.c.j 
+^5abdidi[E{x- - y-)5\x^ - y^)S] , 

= ifahcV^{x)5{x- - y-)S'^{x± - y±) 

-jiifabc + dabc){dt[e{x- - y-)5\x^ - yA_){y^{x\y) - ie'^ A^Mv))] 
+d:,j[e{x- - y-)S\xi_ - yi_)i-SW^ix\y) + ie'^ A+ ix\y))]} ^^'^^ 

- y~)S^ix± - y±)'^{x)j+j'Aabipiy) - /i.c.j 
+^Sabdtdidi[e{x- - y-)S\x^ - y^)S] , 



where V^'* is the usual SU{3) vector current: 

V,^{x)=^{x)rjXaHx), (3.10) 

and V^{x\y) and Ai^{x\y) are the bilocal generahzation of the vector and axial vector cur- 
rents: 

Vr{x\y)=iiix)j^^\ai^iy), (3.11) 

A^{x\y) = ll]{x)j^j^^Xa^l^{y) , (3.12) 

with the internal symmetry matrices that satisfy the relation A^Afe = {ifabc + dabcj^c- 

The fixed mass sum rules based on light-front current algebra are now derived from 
the tensor 

CZiP^ 9) = / d'xe'^^ps\[V,^{x), vamps) . (3.13) 
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The time reversal invariance, Lorentz invariance, parity conservation and current conserva- 
tion allow us to rewrite 



By setting g+ = and integrating over q~ for the +i/ components of C^^ [31], we have 



q+=0 

q~=(i^+pj_-qj_)/p+ 



(3.15) 

= J dx-d^x^e-'P^--Hpsm{x),V,-{0)U=o\ps), 
from which one can immediately obtain the following fixed mass sum rules: 

r duWt^^ {q\ u) = TTfabcrc , (3.16) 

Jo 

/ duW^f\q\ u) = -TTfabc / doAliO, a) , (3.17) 
Jo I Jo 

dyWf'\q\u) = -'Kfabcj^ daAl{<^,a), (3.18) 

r duW'£'\q\v)=Q, (3.19) 
Jo 

r duWt'''\q\iy) =0, (3.20) 
Jo 

roo . . 1 POO , 

/ du{u/-q^)wt'''\q^u) = -7rdabc daF,(0,a), q^ < , (3.21) 

Jo ^ Jo 

where wj^""^^ and W^/"^^ are the symmetric and antisymmetric decomposition of W^''' in ab: 
Wf = W'l"^^ + iWl''^\ and % and are defined as follows: 

(p|A(^(x|0)-A^^(0|x)|p) 

= 2i{si^A^^{x^,x -p) +p'^x ■ sAI_{x'^,x ■ p) + x'^a; • sA^(x^,x ■ p)) , 

(3 22) 

{p\vr{m-vmx)\p) 



= 2i{p'^vl{x'', x-p) + xt'Vlix^ x-p)). 

Eq. (3.16) is the Dashen-Gell-Mann and Fubini sum rule, and Eq. (3.17) the so-called Beg 
sum rule. The sum rules of Eqs. (3.16-21) provide a large number of relations between the 
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various hadronic coupling constants and form factors that can be measured experimentally 
[31]. There are also many other sum rules among the matrix elements of various currents 
that can be derived in the same way. These are the earliest applications of light-front 
dynamics to the strong interaction phenomena. Theoretically, to justify these sum rules, 
one needs to develop the fundamental theory (QCD) on the light-front and to find the 
light-front hadronic states, which allow us to directly compute these matrix elements from 
first principles. 

III-2. Light-Front Axial Charge and Its Dynamics 

In this subsection, I will discuss the properties of the axial current and axial charge 
matrix elements in hadronic states, which are the most important parts for the understand- 
ing of hadronic structure. 

Consider the axial current A^. The ordinary axial charge and the light- front axial 
charge Q^j^ are defined respectively as follows: 

Ql{t) = J d'xA'^ix), (3.23) 

QlLix"-) = \j dx-d'x^Atix) . (3.24) 

Thus, the matrix elements of and Q^j^ in two hadronic states, \a) and |/?), are related 
to the matrix elements of the axial current as follows: 

{l3\Ql{t)\a) = C^^fS^pp - PamAl{t,0)\a) , (3.25) 

{P\Qli{x+)\a) = i27rf S{p+-pt)S\pp^ - p„x)(/3|A+(x+, 0)|«) . (3.26) 

The low-energy hadronic phenomena ensure that the axial current is partially conserved 
(PCAC): 

df,A>i = UmliTa , (3.27) 

where /^r is the pseudoscalar (here I simply call it the pionic) decay constant, ttt-tt the pionic 
mass and vr^ the pionic interpolating field with a = 1, ■ ■ ■ 8 for SU{3) flavor. PCAC implies 
that the matrix element of the axial current is dominated by the pionic mass pole [92,93]: 

{Pm»)m^ = -iUPT^ala) „2 "^^a + (/?l^^l«)iv, (3-28) 



where {p\Ai-^\a) denotes the nonpole matrix element. Using the PCAC and the identity, 
{P\Tra\a) = - 1^-2 {f3TTa\a), we have 

{f3ma) = - fj'^ , )(/3|A-(t,0)|a)^, q = P0-Pa, (3.29) 
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which leads to 

(/3|Q5(0)|a) = {27rf6Hpp - Pa)-^^(/3|AO(0)|a)^ , (3.30) 

{mlim<^) = {2Trf5{p+-pt)5\pp± - ipa±mAtm<^)N ■ (3-31) 

Obviously, the matrix element of in two hadronic states describes the pionic 
transition processes. But in the instant form the matrix element is singular and is zero 
in the chiral limit m^r 0; therefore it does not provide us with much information about 
dynamical chiral symmetry breaking while on the light-front the matrix element is well 
defined and is independent of the chiral limit. Furthermore, the equation of motion for the 
axial charge is 

^^'^ =-[QIl,Hlf]=D, (3.32) 



where H^f = P~ is the QCD light-front Hamiltonian, and D is determined by the above 
commutation. Since the light- front momentum P+, are conserved quantities and com- 
mute with Q\i, this leads to 

[QIl, P+P- - Pl] = [Q, M2] = iP+D , (3.33) 

where is a mass-square operator. Sandwiching the above equation between two hadronic 
states, we have 

{ml - mlmQlLW) = ip-^imo^) ■ (3-34) 

Combining Eqs. (3.31) and (3.34), we have made a simple but direct connection 
between the measurable hadronic phenomena and the fundamental QCD interaction. Obvi- 
ously, the basic theoretical aspect that we need to explore here is H^p and the corresponding 
bound states. 



III-3. Parton Distribution Functions and Fragmentation Functions 

Probing the hadronic structure functions via deep inelastic scattering experiments has 
dominated the study of the strong interaction and hadronic physics over the last twenty-five 
years since the discovery of the scaling rule of the structure functions by Bjorken [12] and 
the interpretation of the scaling behavior via the parton picture introduced by Feynman 
[10] in 1969. The parton picture emerges most naturally on the light-front. In the last 
subsection, we have discussed the sum rules for the hadronic structure functions which 
one extracts from cross sections measured in deep inelastic scattering. Theoretically, using 
the operator product expansion, the hadronic structure functions can be separated into 
hard and soft parts near the light-cone. The hard part is now known as the so-called 
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hard scattering coefficient in QCD, which is the short distance contribution of quarks and 
gluons to the structure functions. The soft part measures the low-energy (nonperturbative) 
properties of quarks and gluons in the parent hadron. When QCD is quantized on the 
light-front with the light-front gauge, the soft contribution can be identified with the parton 
distribution functions, i.e., the number density of partons as a function of the fraction of 
the light-front longitudinal momentum of the parent hadron. In the last two decades, QCD 
has been used extensively to compute the hard scattering coefficients that arc relevant for 
the scale evolution of hadronic structure functions. However, the QCD based exploration 
of the nonperturbative part, the parton distribution functions themselves, is still in a very 
preliminary stage. The current high-energy experiments, such as the e^p HARE experiment 
and pp collisions at FNL and CERN, will provide more accurate information about the 
parton distribution functions over a broader range of the momentum transfer in these 
processes. Currently, the most interesting problem in the study of parton distribution 
functions is to look for the origin of proton spin and the small- a; physics, as well as various 
inclusive heavy meson decay properties. Meanwhile, the study of the hadron productions in 
high energy processes, namely the fragmentation functions which measure the probability 
of finding a hadron in a high energy parton, is also one of the currently most interesting 
subjects. In this subsection, I will discuss the formulation of these problems on the light- 
front. 

1. Distribution functions. In deep inelastic scattering, the cross section is proportional 
to the hadronic tensor which is defined by 

W'^''ip,q) = d^xe'P-{ps\[J''{x),rm\ps) . (3.35) 

We can rewrite W^'^ in terms of p, q, s and the invariant tensors g^'^ and s'^'^'^l^ under the 
constraints of parity, time-reversal invariance and current conservation: 

,) = + i^) + ^ (p. _ ^<,-.) (p^ - J^,^) ^^ ^^^ 

+iW3e^^'''^'^qaSf3 + iW^e^'^^f^ q^iq ■ psp — s ■ qpp) , 

where the Wi's are Lorentz invariant structure functions which determine the underlying 
hadrons. Based on the factorization theorem, the hadronic tensor can be separated into 
hard and soft parts: 

W^^{p,q) = E(-^)' / 4^.}Cr(xB,a(g2),{x,})/,,({x,},-g-) , (3.37) 

where C-t^f is the hard parton scattering coefficient which can be calculated from pertur- 
bative QCD, while fir is defined as the parton correlation function in a hadron, which is 
associated with the nonperturbative aspect of QCD. The notation r represents the twist 
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index, Xi denotes the momentum fraction carried by the i-th parton, is the Bjorken 
scale, the QCD running couphng constant, the momentum transfer by the pho- 

ton, and A the scale of the strong interaction. On the light-front, fir is reduced to the 
parton distribution function: 



where = 1, 7^, 7x7^, nA = x" which means that the two field operators are separated in 
the longitudinal direction. The study of parton dynamics, which has dominated the inves- 
tigation of the strong interaction in both experiments and QCD in the past twenty years, is 
indeed based on the above picture of factorization and parton distribution functions. There 
are a huge number of papers in the literature devoted to this aspect. A nice discussion and 
derivation may be found in Ellis, Furmanski and Petronzio's work [54], and a review paper 
given by Collins, Soper and Sterman could also be a good reference on this topic [55]. Jaffe 




(3.38) 
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TABLE I. 

TwiKt-2: 0(1) Twist-3: 0{l/Q) TwiKt-4: Oil/Q'^) 

Light-front projection 

Target spin average 
Target helicity asymmetry 
Target helicity flip 



1 

7±7^ 



++ 

h{x) 

hiix) 



+- 

e{x) 
h2{x) 
9t(x) 



h{x) 
h^ix) 



and Ji have recently classified various parton distribution functions based on the fermionic 
light-front up and down components and helicity on the light-front. These parton distribu- 
tion functions can be naturally classified in terms of their helicity (spin) dependence and 
twists on the light-front, as we can see from Table L For the detailed definition see Ref. 
[94]. 

2. Fragmentation functions. Fragmentation functions were introduced by Feynmen 
to describe hadron production from the underlying hard parton processes [11]. Collins and 
Soper showed later [95] that, in QCD, parton fragmentation functions can be defined as 
matrix elements of quark and gluon field operators at light-front separations. Physically, 
the quark and gluon fragmentation functions are the probability to find a hadron in a hard 
parton in hard processes, such as the lepton-hadron and hadron-hadron colliders. Explicitly, 
the quark fragmentation is defined by 

fA/i{z) = dx-e-f^--/^TV(O|0+(O)|lf(p)s)(lf(p).|0U^")|O) , (3-39) 

while the gluon fragmentation is 

^""'^'^ = 3^ / dx-e-'^^^- {0\F+>^mH{p)s){H{p)s\F+{x-m , (3.40) 

where Tr traces the color and Dirac components of quarks, \H{p)s) is the hadronic states 
with the momentum p and the helicity s on the light-front. Here we have also taken the 
light-front gauge = 0, F^^^ is the usual gauge field strength tensor. 

3. Inclusive decay of heavy-mesons. Recently it was found that, using the heavy 
quark symmetry and operator product expansion, heavy meson decay spectra can also 
be reformulated in terms of parton distribution functions [97]. Typical examples are the 
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radiative and scmilcptonic i?-mcson decay spectra which are given by 




dT{B XJV) 



^'^parton 



dEi 



{Ei,ml){f{k+) + 0{A/m,)), 



(3.41) 



where — ^^E°" hard 6-parton decay spectra. The universal parton distribution func- 

tion inside the B meson, f{k~^), is defined on the hght-front: 



with /i+ is the hght-front up component of the heavy quark field and \B{v)) is the light-front 
S-meson states [96]. 

In conclusion, we can see that light-front dynamics has played a very important role 
in the study of the strong interaction phenomena in the past thirty years, especially for 
parton physics. Based on the kinematics, symmetry and unitarity of physical strong inter- 
action processes, most of the measurable hadronic structural properties have been reduced 
to matrix elements of light-front field operators in the light-front hadronic bound states 
plus perturbativc QCD corrections. The latter are calculatable from the fundamental the- 
ory. While these nonperturbative hadronic matrix elements (the form factors, the coupling 
constants, the parton distribution functions and parton fragmentation functions) have been 
understood quite well from the experimental data and some phenomenological model stud- 
ies, unfortunately, we still lack a reliable approach to calculate these physical quantities 
from QCD. Obviously, in order to completely understand the strong interaction phenom- 
ena from the fundamental theory, it is natural to require the knowledge of the light-front 
QCD bound states for hadrons. In the next section, I will first discuss the structure of 
light-front wave functions for hadronic bound states. Then, in the subsequent sections, we 
shall look for a possible approach to determine these hadronic bound states from QCD on 
the light-front. 



In terms of the standard language of field theory, relativistic bound states and res- 
onances are identified by the occurrence of poles in Green functions. Although the in- 
formation extracted from this approach provides a definition of physical particles, it does 
not provide their detailed properties, namely their wave functions which determine all the 
structures of the particles. In order to understand the dynamics of hadrons, we need to 
find the explicit form of the hadronic wave functions. 




(3.42) 



IV. STRUCTURE OF LIGHT-FRONT BOUND STATES 
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However, the wave functions of bound states in field theory have not been well 
established. One may define the wave functions as the eigenstates of and may determine 
them by solving the eigenequation of P°. But P° is a square root function of the momentum 
and mass which does not give us a clear picture of the Schrodinger eigenstate equation in 
quantum mechanics. The well-defined framework for finding the relativistic bound states is 
to solve the Bethe-Salpeter equation. However, the Bcthc-Salpctcr equation only provides 
the amplitude of a Fock sector in the bound states so that it cannot be normalized. In other 
words, the Bethe-Salpeter amplitudes do not have the precise meaning of wave functions 
for particles. 

In principle, the wave function of a relativistic bound state can always be written as 
an operator function of the particle creation operators acting on the vacuum of the theory. 
However, for many theories that we are interested in, especially for QCD, the vacuum is 
very complicated. Thus, it is even more difficult to formally write down a relativistic bound 
state wave function. 

In this section, I will discuss how these difficulties may be solved when we look at 
the problem on the light-front. 

IV- 1. Light-Front Vacuum 

In the equal-time framework, the vacuum of QCD is crucial for a realization of 
chiral symmetry breaking. From the vacuum the axial charges (a = 1, • • • 8) can create 
pseudoscalar particles which are the lowest bound states in the strong interaction regime 
that have been used as the building blocks in the development of the effective chiral field 
theory for low-energy hadron dynamics. Undoubtedly, the complicated vacuum is the most 
important object to be understood for the nonperturbative aspect of QCD. As we have 
already argued, it is also a starting point for the construction of hadronic wave functions. 
However, the understanding of the true QCD vacuum is still very limited, although a lot 
of informative work has been carried out in the past two decades based on the instanton 
phenomena [66] and the QCD sum rule [71]. 

In the light-front coordinates, a particle's momentum is divided into the longitudinal 
component (along the light-front time direction = t + x"^) and the transverse components 
(perpendicular to the light-front time direction). The longitudinal momentum of each 
particle, p+ = + p^, cannot be negative since the energy of a physical state always 
dominates its momentum. As a result, the light-front vacuum for any interacting field 
theory can only be occupied by the particles with zero-longitudinal momentum, namely 

\vac)LF = f{al+^oM, (4.1) 

so that P~^\vac)LF = 0, where P+ = Ylii^f ■ At this point, the light-front vacuum is 
still nontrivial. In the past several years, many tried to solve the so-called zero-mode (the 
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particles with = 0) problem to obtain a nontrivial light-front vacuum but not much 
useful progress has been made since then [98-107]. 

To construct hadronic bound states in the Fock space of quarks and gluons, it is 
natural to ask whether one can simplify the vacuum so that we can easily construct the 
hadronic states by the use of the Fock space expansion with a trivial vacuum. It is obvious 
that if we can remove the basic constituents with zero longitudinal light-front momentum, 
the vacuum of the full interacting theory is the same as the free field theory, namely 



To remove these particles with zero longitudinal momentum, we can either use a 
prescription that requires the field variables to satisfy the antisymmetric boundary condition 
in the light-front longitudinal direction [88] or deal with a cutoff theory that imposes p+ > e 
on the momentum expansion of each field variable, where e is a small number [85]. Thus, the 

positivity of longitudinal momentum with such a prescription or an explicit cutoff ensures 
that the light-front vacuum must be trivial, as has been shown above. Now the wave 
function of a relativistic bound state can be expressed as an ordinary Fock state expansion: 



longitudinal momentum, and /(a^) must be a color singlet operator. 

Of course, in this case, the light-front axial charge, denoted by Q^l, vanishes the 
vacuum state. One might ask that, once the vacuum is simple, what happens to the 
nontrivial effect of the theory associated with the nontrivial vacuum. At least we believe 
that confinement and chiral symmetry breaking in QCD must be a result of the nontrivial 
vacuum. In the next subsections, we will take two simple examples, the sigma model and 
1-1-1 QCD, to show how the zero-mode particles are removed but the nontrivial vacuum 
effects and confinement mechanism can be recovered or manifested. But before we explore 
the nontrivial vacuum effects with the trivial vacuum structure, I would like to discuss first 
the light-front bound state equation. 

IV-2. Light-Front Bound State Equation 

Once the light-front vacuum becomes trivial, the light-front bound states for various 
hadrons can be expanded in terms of the Fock space. Explicitly, a hadronic bound state 
labeled by a with the total longitudinal and transverse momenta P+ and P±, and the 
helicity (the total spin along the longitudinal direction) A can be expressed as follows: 



vac)LF = jO) . 



(4.2) 




dxi(Pk±i 



n, XiP+,XiP± -h k±i, Xi)^j^/a{xi, k±i, Aj) , (4.4) 



a 



2(27r)3 
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where n represents n constituents contained in tire state |n, XiP^ , XiP± + k±i, Aj), Aj is the 
hehcity of the i-th constituent, and /' denotes the integral over the space: 

^ Xj = 1 , and XI ^-Lj = ^ ' (^■^) 

i i 

here Xi is the fraction of the total longitudinal momentum that the i-th constituent carries, 
and k_\_i is its relative transverse momentum with respect to the center of mass frame: 

Xi = , ki±_ = pi± - XiP± , (4.6) 

with pf, pi± being the transverse and longitudinal momentum of the i-th constituent, 
and ^n/ai^i, k±i, Xi) the amplitude of the Fock state \n,XiP^ ,XiP± + k±i,Xi) with the 
normalization condition: 

E f ^^S^\^n/a{XiMiAi)f = 1. (4.7) 

The eigenstate equation that the wave functions obey on the light-front is obtained 
from the operator Einstein equation P^ = P~^P~ — = M^: 

HLF\a, P+,P±, A) = -L^^\a, P+,P±,X), (4.8) 

where Hlf = P~ is the light-front Hamiltonian. Futhermore, since the boost on the 
light-front only depends on kinematics, we can consider the bound state in the rest frame 
{P± = 0, P+ = Ma). Thus, the eigenstate equation simply becomes: 

HLF\a,P^,P±=0,X) = Ma\a,P+,P±=0,X), (4.9) 

which is the familiar Schrodinger equation in ordinary quantum mechanics. On the light- 
front, boosting a bound state in the rest frame to any other frame is quite simple, and 
is dynamically independent, as we have shown in Eq. (2.10). Thus, once we find the 
bound state in the rest frame, we can completely understand the particle structure in any 
frame. Yet this is not true in the instant form. In the instant form, although the bound 
state equation in the rest frame has the same form, the solutions in the rest frame are 
not easily boosted to other Lorcntz frames due to the dynamical dependence of the boost 
transformation. Therefore, in each different Lorentz frame, one needs to solve the bound 
state equation of P^ to obtain the corresponding wave functions. This is perhaps the reason 
why we have not established a reliable approach to construct relativistic wave functions in 
the instant field theory in terms of the Schrodinger picture. This obstacle is obviously 
removed on the light-front. 
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To sec the explicit form of the hght-front bound state equation, let us consider a 
meson wave function (for instance, a pion, [47]). The light-front bound state equation can 
be expressed as: 



m 



2 V- +mf 



11 



qqg 



/ {qq\Hint\qq) {qq\Hint\qqg) 
{qqg\Hint\qq) 

V 



(4.10) 



Of course, to exactly solve the above equation for the whole Fock space is still impossible. 
Currently, two approaches have been developed. One is given by Brodsky and Pauli [108- 
111], the so-called discretize light-front approach, the other by Perry, Harindranath and 
Wilson [112-115], based on the old idea of the Tamm-Dancoff approach [116,117] that 
truncates the Fock space to only include these Fock states with a small number of particles. 
Furthermore, if one can eliminate all the high order Fock space sectors (approximately) by 
an effective two-body interaction kernel, the light-front bound state equation is reduced to 
the light-front Bethe-Salpeter equation: 



kj_ + 
x(l ~ x\ 



^qq{x,k±) = j 



dycfk^ 
2(27r)3 



Veff{x,k±,y,k'_^)-^qq{y,k'_^) . 



(4.11) 



The advantage for Eq. (4.10) with the Tamm-Dancoff approximation is that it provides a 
reliable way to study the contribution of Fock states which contain more particles step by 
step by increasing the size of truncated Fock space, while the Bethe-Salpeter equation lacks 
such an ability. 



IV-3. Spontaneous Symmetry Breaking on the Light-Front 

Now, we turn back to the question that if the vacuum on the light-front becomes 
trivial, how are we to realize the physics of the spontaneous chiral symmetry breaking 
and the associated Goldstone bosons which are all the theoretical underpinning for our 
understanding of the low-energy theorem. We argue that eliminating the longitudinal zero 
modes leads to a trivial light-front vacuum, as a result, the theory must become explicitly 
symmetry breaking by a proper modification of the Hamiltonian. To be explicit, let us 
consider the sigma model as an example. The sigma model is one of the few models in field 
theory that can clearly demonstrate the mechanism of spontaneously symmetry breaking. 

The Lagrangian density of the sigma model is given by 

C = ^{{di^af + {d'^Trf) - Via'' + tt^) , (4.12) 
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where 

V{a^ + TT^) = + + + ■ (4-13) 

It is well-known that this Lagrangian has a continuous rotational symmetry in the u — tt 
plane. This symmetry is spontaneously broken since the vacuum expectation value of the 
field variables, for example the a field, does not vanish (the vacuum is nontrivial). As a 
result the vr is a Goldstone boson. 

Now we want to define the theory on the light-front such that the vacuum is trivial 
but the massless tt can be determined in an alternative way (i.e., it is not realized from the 
mechanism of spontaneous symmetry breaking) . 

Since the original theory contains a zero mode in the a field, we can remove the zero 
modes by shifting the a field to a' such that cr' does not contain zero modes: a' = a — a^. 
The resulting potential in terms of ((T',7r) is 



V[a' ,Tx) = \[—{a 4- TT ) + (7„cj (cr + vr ) + — a^a + —a^-n + a^a + 



(4.14) 



With this potential, the vacuum is trivial. Now the question is, what is the value of cr^ 
which cannot be determined in the theory with the zero modes removed. However, current 
conservation can be used to fix a,^J. The continuous transformation of the original potential 
leads to a conserved current: = ad^'K — ird'^a. After the shift, it becomes: 

f = a'd^^TT - ndi^a' + a^d^'ir . (4.15) 

The dropping of the zero longitudinal momentum modes however does not affect local 
current conservation, d^j^ = 0. It is easy to check that current conservation results in the 

/ 2 

value of ay = J 4r . The potential then is reduced to 



V{a',n) = \f.^a^ + ^ia'^ + .r + ]l^ (a ^ + tt V , (4.16) 

which shows that the pion is massless. One can now find that the charge Q = ^ J dx~(Px±j~^ 
no longer commutes with the Hamiltonian. Thus, in the theory with zero modes dropped, 
the Hamiltonian explicitly breaks the symmetry and the vacuum is trivial but the pion is 
still massless. 

It may be noted that the above analysis is also true for the instant formulation but 
the conclusion is different. In the instant form, we can follow the same route (by a shift of 
the a field) to determine the symmetry breaking and the massless vr but it has nothing to 
do with the simplification of the vacuum. The shift in the instant form does remove the 
zero modes from the theory but the instant vacuum is not only occupied by the zero modes. 
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Of course, the above demonstration of symmetry breaking, with a trivial vacuum 
on the hght-front, cannot be directly applied to QCD. The situation with QCD on the 
light-front is much worse than the sigma model. In QCD, the operator that is expected to 
have a non-zero vacuum expectation value due to spontaneous symmetry breaking is the 
composite operator i/jtp where %p is the quark field. There is no obvious mechanism for this 
expectation value to show up in the QCD Hamiltonian since the QCD Hamiltonian does not 
apparently contain squares or higher power terms of ^p■^JJ. This is perhaps the reason that all 
these beautiful but naive pictures of spontaneously symmetry breaking are not manifest in 
QCD. However, the trivial vacuum with an explicitly symmetry breaking Hamiltonian on 
the light-front may provide a new avenue to explore the physics associated with Goldstone 
particles. Therefore, for QCD, the light-front zero modes are not the essential problem we 
need to solve. The problem is how we can obtain a correct effective Hamiltonian after the 
elimination of zero modes. 



IV-4. Bound States in 1+1 QCD and Confinement 

Perhaps a clear understanding of quark confinement is best seen in the quark dynamics 
of 1+1 QCD on the light-front [118]. Here we will demonstrate the confinement mechanism 
on the light-front in this theory with a trivial vacuum. The Lagrangian density for 1+1 
QCD is given by 

£1+1 = -iFr^aM. + ?(i7 • 1? - m)V , (4.17) 

where F^^" = d^A"^ - d^Ai^ - gP^'^A'^A'^ and = + igT^'A'^, with /x, 1/ = + or -. If 
we take the light-front gauge, A'^ = 0, the Lagrangian can be rewritten as 

A+i = i^id-^ + + gjtK - \{d+A-f , (4.18) 

where V'- and A~ are the light-front constraint variables that satisfy the following constraint 
equations: 

= -im(3-^ij+ , (4.19) 

d+A- = 2g^j+, j+ = 2V'tr>+. (4.20) 
Then the light-front 1-|-1 QCD Hamiltonian is simply given by 

/^i+i = J dx-{m2^^^V+ - 2g'jt-^jt} ■ (4.21) 

The confinement can be seen explicitly if we use the definition of Eq. (2.16) for the operator 
In such a definition, we have removed the zero longitudinal modes. The Hamiltonian 
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becomes 



i^i+i = — J dx~^dx' i^—im'^ip\.{x )e{x — x' )tl)^{x' ) 

-g^j+ix-)\x--x'-\j+{x'-)] + Qim A^(y dx-j+Cx")))' 



(4.22) 



In this Hamiltonian the first term is the free quark mass term, the second term is a hnear 
confinement potential between quarks, and the last term is a singular (divergent) term due 
to the elimination of the zero modes. It is this term that forbids the observation of a single 
quark. To see how the confinement mechanism is manifested in this theory, let us consider 
the color non-singlet single quark state and the color-singlet two quark bound state. 

Let us first consider a single quark state: |P+,c) = 6j(P+)|0), where 5j(p+) is the 
quark creation operator with momentum j)+ and color index c; 

V'c+(x-) = I ^[6e(p+)e^f^-"/2 + 4(p+)e-^f^-"/2] , (4.23) 

satisfies the light-front fermion commutation relations. It is easy to check that the solution 
of the single quark bound state equation, 

P-|P+,C) = H,+,\P+,c) = -^\P+,c) , (4.24) 



IS 



o2 2 p+ 

with Cf = T^-T" = |, e = ^ ^ 0(A ^ oo). Obviously ^ oo (e ^ 0) which 
implies that no single quark can be measured. In other words, the quark is confined. This 
confinement mechanism is obviously an effect of the original nontrivial vacuum since the 
term ^ originates from the elimination of the /c^ = singularity. 
Now we consider the quark-antiquark color singlet state: 

= / ^$(p+)6t(p+)4(P+ - p+)|0) . (4.26) 
The eigenstate equation of the above bound state is 

7 r = —Cf / -y- J-'- . (4.27) 



x(l — x) J TT J At: {x — y)"^ 

Here the confinement of two quarks is given by a linear potential (the right hand side of 
the above equation). The singular term in Eq. (4.22) does not contribute to the above 
bound state equation since it is proportional to the color charge operator, and should 
be finite. 
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The above equation does not have an analytic solution. The discussion about the 
general form of the wave function can be found in Ref. [118]. The approach can be 
applied to 1+1 QED where an additional property, namely the gauge anomaly, can be 
addressed canonically [119]. However, I must emphasize that all the conclusions about 
quark confinement obtained here disappear when we consider the real theory of the strong 
interaction, namely 3-f-l QCD. As we will see later, in 3-1-1 QCD there are physical gauge 
degrees of freedom (the transverse gaiige field) which cancel the effect of confinement from 
the linear potential due to the elimination of the unphysical gauge degrees of freedom. 
Nevertheless, here it shows that eliminating the zero longitudinal momentum modes leads 
to a trivial vacuum but the confinement mechanism can be manifested in the Hamiltonian. 
Undoubtedly, the manifestation of confinement in 3+1 QCD with a trivial vacuum must be 
much more complicated and we will address it later. 

It may be worth mentioning that in the last few years, many studies on nonperturba- 
tive features of light-front dynamics were focused on the 1+1 field theory. Typical examples 
are: the discretized light-front quantization approach for the bound states in the 1 + 1 field 
theory developed by Pauli and Brodsky [108,120], spontaneous symmetry breaking and vac- 
uum structure in the 1+1 Yukawa model first considered by Harindranath and Vary [121], 
the light-front Tamm-Dancoff approach for bound state Fock space truncation discussed by 
Perry et al. [112,113], and other extensive discussions of various 1+1 field theory models 
of these aspects [122-128]. 

IV-5. Phenomenological Hadronic Bound States on the Light-Front 

At the present time, how to solve for the bound states, discussed above in this section, 
from 3+1 QCD is still unknown. Hence, it may be useful to have some insights into the 
light-front behavior of the meson and baryon wave functions which have been constructed 
phenomenologically in describing hadrons. In fact, the phenomenological light-front meson 

and baryon bound states have been studied extensively in the last few years, this is called 
the relativistic quark model which was initially proposed by Terent'ev in 1976 [129,130]. 
The motivation is to provide a simple relativistic constituent quark model for mesons and 
baryons that can yield a consistent description of the hadronic properties for both low and 
high Q^. 

The general construction of the phenomenological wave functions is based on the 
non-relativistic constituent quark model. The constituent quark model has been very suc- 
cessful in the description of hadronic structure based on a very simple structure, namely 
that all mesons consist of a quark and antiquark pair and the baryons are made of three 
constituent quarks, their wave functions satisfy the SU (6) classification and Zweig's rule 
which suppresses particle production in favor of rearrangement of constituents for hadrons 
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[76-78,157,80]. However, such a simple picture is very difficult to understand from QCD, 
due to its nonrelativistic assumption and due to our belief that the QCD vacuum must be 
very complicated so that hadrons must contain an infinite number of quark-antiquark pairs 
and gluons. 

The light-front wave functions describe the relativistic hadronic structure with a 
nonrelativistic form. Furthermore, the simple vacuum on the light-front ensures the validity 
of the Fock state expansion of hadronic wave functions. With the assumption of existence 
of constituent quarks (with masses of hundreds of McVs), the leading approximation to 
hadronic wave functions that consist of a quark-antiquark pair for mesons and a three- 
quark cluster for baryons should be a reasonable starting point. 

However, it is not so easy to identify the light-front hadronic wave functions with 
hadronic states. The difficulty is their spin structure. On the light-front, we are unable to 
kinematically construct the hadronic wave functions with fixed spin. The light-front wave 
functions discussed in the last section are labeled by helicity not spin. In these calculations 
of the parton distribution and fragmentation functions, the hadronic bound states arc de- 
fined or classified in terms of the helicity (=chirality). However, when we use the light-front 
wave functions to compute the hadronic structural quantities, such as hadronic decay form 
factors and coupling constants, we must have states with a definite spin. A general solution 
to the spin problem on the light-front has not been found. However, phenomenologically, 
the helicity part of the wave functions on the light-front can be transformed to a light-front 
spin wave function part via the so-called Melosh transformation (which is exact only for 
free quark theory) such that the hadronic states can be projected (approximately) from the 
set of light-front wave functions with fixed helicity [131,132]. Here, I list the detailed meson 
and baryon light-front wave functions that have been used to calculate various hadronic 
quantities in the past few years. 

The general form of the phenomenological light-front hadronic bound states has a 
similar structure to the constituent quark model states: for meson states (Eq. (4.4) with 
only the qq Fock space sector). 



\P+,P±,SS3) 




J2 '^^%x,k±,Xi,X2)\x,k±,Xi;l-x,-k±,X2) 



(4.28) 



and for baryon states (Eq. (4.4) with the three quark Fock space sector). 




X \xi,ki±,Xi;x2,k2±,X2; I - xi - X2, -{ki± + k2±), A3) , 
where ^^^^ is the amplitude of the corresponding qq or three quark sector (the wave function 
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of the quark model): 

^553 ^ jr^SSs^j^,^^ Xi)^{xi, ki±) , (4.30) 

with JT the flavor part of the wave function which is the same as in the constituent quark 
model, and E and $ are the spin and space parts that depend on the dynamics. By using 
the Melosh transformation, 

rrii + XiMo - ia ■ {n X ki±) 
RM{ki±,mi) = , — , (4.31) 

where n = (0, 0, 1), a is the Pauli spin matrix, rrii the i-th constituent quark mass, and Mq 
satisfies 

Xi 



M'o = E ^^^^ ' (4-32) 



the light-front spin wave function is then given by 
Ei^^{k^,Xi,\2) 

= E {>'i\RMik±,mi)\si){X2\RMi-k±,m2)\s2)l^si, ^S2 



SSf 



(4.33) 



for mesons; for baryons the spin part is rather complicated for a detailed construction, see 
Ref. [133]. The momentum part of the wave function is 

^m{xi, h^) = Mm eM-Mi /2pl) , (4.34) 

for mesons and for baryons 

M^i, ki±) = Mb (i^j^^2/^2)3.5 ' (4-35) 

where M is the normalization constant and /3 is a parameter fixed by the data. The 
exponential function for the momentum wave function of baryons has also been used but 
Schlumpf found that the pole-type wave function is much better fit to the data with high 
momentum transfer [133]. 

These phenomenological light-front wave functions have been widely used to calculate 
hadronic form factors and coupling constants; the results are in pretty good agreement with 
experiments for a very broad range of momentum transfer; (see Fig. 3 and 4), and are 
much better than the nonrelativistic constituent quark model and other phenomenological 
descriptions [134-137,133]. 

There is also another phenomenological light-front wave function derived by Chernyak 
from QCD sum rules [138]. Such a wave function is valid for relatively high momen- 
tum transfer, and recently it has also been applied to describe various hadronic processes 
[97,139,140]. 
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However, all these arc just a phenomenological understanding of light-front hadronic 
wave functions. If these wave functions are very close to the real wave functions of hadrons, 
they must be solved dynamically from QCD. In the remaining sections of this article, I will 
discuss QCD on the light-front. 
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FIG. 3. Proton form factor calculated by the use of the phenomenological light-front proton wave 
function [133] . The solid line is given by the pole-type wave function while the dashed line 
is based on the Gaussian-type wave function. 

V. LIGHT-FRONT QCD 

As we have mentioned several times, an intuitive physical picture of high-energy pro- 
cesses is provided by the partonic interpretation [11]. We have already seen that such a 
picture emerges most naturally in the light-front canonical quantization with the light-front 
gauge A'^ = + = 0. In the last two decades, QCD with either the covariant gauge or 
the light-front gauge has been used extensively to compute the hard scattering coefficients 
that are relevant for the scale evolution of hadronic structure functions. However, the QCD 
based exploration of the nonperturbative part, such as the parton distribution functions, 
is still in a very preliminary stage. One may also note that many practical calculations for 
short distance QCD have been performed in the standard Feynman perturbation theory 
with the light-front gauge [42,141,142,54,143]. In this scheme, Feynman rules are derived 
by the use of the path integral approach with a (light-front) gauge-fixing term [40] and var- 
ious hard scattering coefficients are calculated via either the operator product expansion 
[42,141] or a Feynman diagrammatic approach [142,54,143]. Despite the successful per- 
turbative calculations, it is still not clear whether one can extend the Feynman approach 
to nonperturbative studies, such as the QCD calculation of parton distribution functions, 
which requires the full information of the hadronic bound states. 

Current attempts to explore nonperturbative QCD on the light-front, as we have em- 
phasized throughout this paper, are based on the Hamiltonian theory, which is defined by 
quantizing the theory on a light-front plane via equal-x"*" (the light-front time) commutation 
relations, as was first developed by Kogut and Soper for QED [21]. In the light-front QCD 
Hamiltonian theory, the hadronic bound states may be obtained by diagonalizing a light- 
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FIG. 4. Pion form factor calculated by the use of the phenomenological light-front pion wave 
function [134]. The different lines correspond to different Parameters used in the phe- 
nomenological wave function. 

front QCD Hamiltonian in a truncated quark-gluon Fock space based on the old ideas of 
Tamm and DancofF [116,117]. In this formalism the low-energy hadronic structures, namely 
the parton distribution functions, parton fragmentation functions, various hadronic form 
factors and strong interaction coupling constants, can be addressed directly from QCD. 
Meanwhile, the short distance behavior can also be studied in the same framework, i.e., 
in the old-fashioned perturbation theory or explicitly in the a; ''"-ordered perturbative QCD 
Hamiltonian theory based on equal-x''" commutation relations [47]. Indeed, the interpre- 
tation of high-energy processes via the parton picture had led to extensive investigations 
of perturbative field theory in the infinite momentum frame in the early 70's, as wc have 
mentioned in the introduction. At that time Drell, Levy and Yan [15] and Bjorken, Kogut 
and Soper [17] had already pointed out that, in the old-fashioned theory, the physical pic- 
ture for various real physical processes becomes much clearer. The parton picture is just a 
typical example. 

Yet, beyond applications to the exclusive processes given by Lepage and Brodsky in 
the early 80's [47], the x"'"-ordered Hamiltonian QCD has not been explored and exten- 
sively utilized in the last decade. Only very recently, loop calculations in the ^'''-ordered 
perturbative QED and QCD theory have been performed [144-147,82,148]. It is seen that 
there are severe divergences in light-front QCD which are associated with light-front gauge 
singularities. However, a systematic computational method, which involves various regular- 
ization and renormalization schemes for light-front singularities and ultraviolet transverse 
divergences, has not been well established. In fact, although the light-front gauge is known 
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to be a convenient gauge in practical QCD calculations for short distance behavior, there 
are persistent concerns about its utility because of its "singular" nature. The study of 
nonperturbative QCD on the light-front for hadronic bound states requires one to gain a 
priori systematic control of such gauge singularities. Clearly, it will be very helpful to first 
understand the canonical structure of light-front QCD and the origin of various singular 
problems before we go on to explore the nonperturbative properties light-front QCD. 

In this section I will discuss the canonical form of light-front QCD, the origin of the 
light-front gauge singularity, the light-front two-component formulation of QCD, the old- 
fashioned perturbation theory, various severe infra-red (IR) divergences occurring in the 
old-fashioned light-front Hamiltonian calculations for QCD and the associated nontrivial 
QCD structures. I will also discuss some of the difficulties caused by the light-front gauge 
singularity in applications to both the old-fashioned perturbative calculations for short dis- 
tance physics and the upcoming nonperturbative investigations for hadronic bound states. 

V-1. Canonical Light-Front QCD 

The QCD Lagrangian is 

jC = -^TriF'^'^Fi,,) +^{iji,D^' - m)^ , (5.1) 

where F'"' = di^A" - d^Ai^ - ig[Ai^,A% Ai" = Ea^a^" is a 3 x 3 gluon field color 

matrix and the T" are the generators of the SU{3) color group: [T",T^] = if°-^'^T'^ and 
Tr(r'*r^) = i^Sab- The field variable describes quarks with three colors and Nf flavors, 
= — igA^ is the symmetric covariant derivative, and m is an Nf x Nf diagonal quark 
mass matrix. The Lagrange equations of motion are well-known: 

d^Fr + gr''A,^Fr + #7"r> = o , (5.2) 

-m + g-f^^A^")^ = . (5.3) 

1. Light-front QCD Hamiltonian. The canonical theory of QCD on the light-front 
is constructed as follows. The conjugate momenta of the field variables {Ai^{x),ip{x)} are 
defined by 

In terms of the field variables and their conjugates, the Lagrangian density can be rewritten 
as follows: 

= [^F+'id-A",) + zV'U^"V'+)} -n- {A-Ca + i^lC} , (5.6) 
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where 



(5.7) 



and 



(5.8) 



C = (i5+ + gA'^)ip_ - {ia± ■ d± + ga± • yl_L + /3m)ip+ . 



(5.9) 



In Eq. (5.7), we have defined B~ = F^^ as the longitudinal component of the hght-front 
color magnetic field. 

The reason for writing the Lagrangian in the above form is to make the Hamiltonian 
density and also the dynamical variables and constraints manifest. Here the time derivatives 
have explicitly been separated from the others, from which it immediately follows that only 
the transverse gauge fields vl^ and the up-component quark fields V'+ and V'+ are dynamical 
variables. The Hamiltonian density H contains three parts; the first part involves the light- 
front color electric and magnetic fields; the second, the usual quark Hamiltonian with 
coupling to the gauge field and the last a surface term. The last terms in this new form of 
the QCD Lagrangian indicate that the longitudinal gauge field A'^ and the down-component 
quark fields are only Lagrange multipliers for the constraints Ca,C = 0. This is true even 
for the free field theory, as we have already seen in Sec. 11. With the use of the equations of 
motion, one can verify that the gauge field constraint, Ca = 0, is in fact the light-front Gauss 
law which is an intrinsic property of gauge theory. The fermion constraint, C = 0, is purely 
a consequence of using the light-front coordinates. The existence of constraint terms simply 
implies that QCD in the light-front coordinates is a generalized Hamiltonian system. These 
constraints are all secondary, first-class constraints in the Dirac procedure of quantization 
[86]. To obtain a canonical formulation of light-front QCD for non-pcrturbative calculations 
we need to explicitly solve the constraints, namely to determine the Lagrange multipliers, to 
all orders of the coupling constant. Generally, it is very difficult to analytically determine 
the Lagrange multipliers from the constraints Ca,C = since they are coupled by A'^. 
But if we choose the light-front gauge, A'^{x) = 0, the constraints are reduced to solvable 
one-dimensional differential equations: 



^d^E- = d'El + gir^'^bEl - V'lT'^V+) 



(5.10) 



id'^^ijj- = {ia± ■ d± + ga± ■ A± + f3m)il)jf. . 



(5.11) 
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Now the light-front QCD Hamiltonian can be expressed in terms of the physical 
degrees of freedom: 

•,2 



H = 

+ 



(5.12) 



+4 (^) (/"''45+4 + 2^lr«^+) 

-|-surfaceterms, 

where the operator ^ is defined by Eq. (2.16), and the surface terms are given by the last 
term in Eq. (5.7). 

2. A phase space approach to light-front quantization. A self-consistent canonical 
formulation requires that the resulting Hamiltonian must generate the correct equations 

of motion for the physical degrees of freedom (j4Jj, ^/j+j ^^). To reproduce the Lagrangian 
equations of motion, we need to find consistent commutators for physical field variables. In 
the light-front gauge, the (symmetricalized) Lagrangian of Eq. (5.6) can be rewritten as 

jr = 1(5+45- + #i5-V+ - i{^-^|;i)^+) - n . (5.13) 

Thus the light-front QCD phase space is spanned by the field variables, A^^, and their 

canonical momenta, = ^5+^1^, 7r^_^ = ^tp\., vr^t = ~^'^+- The phase space structure 
which determines the Poisson brackets of its variables can be found by the Lagrangian 
one-form Cdx'^ (apart from a total light-front time derivative) , 

Cdx+ = ^2{SidAi + Tr^^dij+ + dVUv,; ~ ^""^^^ ~ t^7rv,+V+ - V^+^^tt^O 

-Hdx+ (5.14) 
= ^q''T^pdqf'-ndx+, 

where the first term on the right-hand side is called the canonical one-form of the phase 
space, and quark fields are anticommuting c- numbers (Grassmann variables). Correspond- 
ingly, the symplectic structure or the Poisson bracket of the phase space is given by 

CO = ^r^pdq-dq^ or [q^, g"]^ = F^^ . (5.15) 
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Canonical quantization is then obtained by replacing the Poisson brackets by the equal-x"*" 
commutation relations: 

[q^{xU"{y)U=y+=iV-l. (5.16) 

Explicitly, we have 

Wa{^)A{y)U=y^ = i^Sa,6^'S\x - y) , (5.17) 

{V'+(x),7r^^(y)},+^^+ = i^5\x - y) , (5.18) 

{^l{x),7r^,^{y)}^^^^^ = -z^S%x-y), (5.19) 
or simply 

[Al{x),Al{y)l^^^^ = -i^e{x- - y-)5\x^ - y^) , (5.20) 

{^+{x),i^i{y)},+^y+ = zA+53(x - y) . (5.21) 

Prom these commutation relations it is straightforward to verify that the Hamiltonian 
equations of motion are consistent with Eqs. (5.2) and (5.3): 

= i^igA- - {a_L • {id± + gA_i_) + Pm} (5.22) 



d-Ai= j[AlH] 



-^[DLFt - Dl,E, - gjl - gr'^A^d+Al] , 



(5.23) 



where D^^ = Sabd' - gP'^'Ai. 

V-2. Light-Front Gauge Singularity 

Note that the above canonical formulation docs not completely define theory for 
practical computations due to the existence of the gauge singularity. The gauge singularity 
is perhaps the most difficult problem in non-abelian gauge theory that has not been com- 
pletely solved since it was developed. In light-front QCD, it arises from the elimination of 
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the unphysical gauge degrees of freedom. To eliminate the unphysical degrees of freedom 
on the hght-front, we need to solve the constraint equations which depend on the definition 
of the operator In Sec. II, we defined this operator by Eq. (2.16). In general, we 

have 

/ 1 \ 1 f°° 

\^'Q+)fi^~^^^'^-^) = ^J dx^£{x~ - x];)f{x];,x+,x±) + C{x+,x±) , (5.24) 

where e{x) has been given in Eq. (2.16), and C{x'^,x±) is a x~ independent constant. The 
reason for choosing the definition of Eq. (2.16) is that the light-front initial value problem 
is associated with the longitudinal boundary condition. A suitable definition of l/c?"*" can 
uniquely determine the initial value problem at x"*" = for independent field variables [32]. 
Rohrlich has shown that such a suitable definition is to set C{x~^,x±) = [i.e, Eq. (2.16)] 
[22]. This corresponds to choosing an antisymmetric boundary condition for field variables 
in the light-front longitudinal direction. I must emphasize here that this is not necessarily 
the best choice for the operator As we will see later this definition still causes many 
undesirable problems associated with the gauge singularity, which even exist in perturbative 
calculations. It is still an open question whether one can find a better definition to solve 
the light-front gauge ambiguity in Hamiltonian formalism. 
By using Eq. (2.16) and the following identity [149], 



— / dx s{x — x )e{x — x" ) = —\x — x" |H A, (5.25) 

2 y-A/2 2 



where the parameter A denotes the distance between two boundary points in the longitudinal 
direction, the constraint gauge and quark components can be determined: 

E-{x) = -d'Alix) - I r dx'-E{x-x'-)jt{x+,x'-,x^) , (5.26) 
Aa{x)= dx'-[e{x- -x'-)d'Al{x+,x'-,x^) 

^ J —oo 



+g\x -x' |j+(a;+,a;' , x_l)] + ( lim A ) -^j+(a; ,x) 

A— >oo / 4 



(5.27) 



il)-{x) = — — / dx d^x'^e{x — x )5'^{x^—x'^) 
4 J— oo 

X [a_L • + gA^{x')) + (5m]'il)+{x') , 



(5.28) 



where j+ = 2'(/'^_T"'^_|_. In this solution, the first surface term at the longitudinal infinity in 
Eq. (5.7) vanishes. Moreover, it is reasonable to assume that the transverse color electric 
fields El^ as well as A^^ vanish as 0(r~^) and 0{r~^) at r = |a;j_| — oo because the gauge 
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freedom is totally fixed at transverse infinity. Thus the second surface term vanishes as 
well. Now the light-front QCD Hamiltonian is simply given by 

= j dx-d\^i^^{d^A^,f + gr'^AlAid^Ai + I- r^"^ f-^^^ A^AiAl^Ai 

+- / dx'- [2gd'A\e{x- - x'-)pa{x'- , x) 

^ ■'-'^ (5.29) 
-iijX{a\id' + gA^) + I5m}e{x- - x'-){a^{id^ + gA>) + /3m}V'+] 

■ / dx'- Pa{x- ,x)\x- - x'-\pa{x'- ,x)\ 
J —oo J 

-l-Qim J d^x±^^J dx-pa{x-,x)^ . 

From Eq. (5.29), we see that, after eliminating the longitudinal gauge field in the light-front 
gauge, a color charge instantaneous interaction emerges in the Hamiltonian, with a linear 
potential in the longitudinal direction. In addition, Eq. (5.29) also contains a singular 
boundary term (the last term) associated with the linear instantaneous interaction. There 
is no such term in the Coulomb gauge because the Coulomb potential vanishes at spatial 
infinity. This term is very similar to that in the 1 -|- 1 QCD result of eliminating 

the unphysical gauge degrees of freedom. However, the role it plays here is totally different 
from that in 1-|-1 QCD. As we will show later, in perturbation theory this term is regularized 
by the distribution function of the product of two principal value prescriptions and leads 
to the cancellation of the light-front linear infrared divergences. For physical states, the 
requirement of the finiteness of energy results in asymptotic equations for the transverse 
gauge fields, which shows that the asymptotic transverse gauge fields do not vanish at 
longitudinal infinity, as we shall sec later. 

The gauge singularity manifests itself more clearly in momentum space. In momentum 
space, the gauge singularity means that the constraints Eqs. (5.8) and (5.9) do not uniquely 
determine the dependent fields in terms of physical fields at longitudinal momentum k'^ = 0. 
In other words, simply choosing A'^ = with the definitions of Eq. (2.16) for ^ does not 
completely fix the gauge degrees of freedom [150]. As we have mentioned many times 
throughout this paper, the k'^ = particles are the constitTients leading to a complicated 
vacuum. From the momentum representation of Eq. (2.16), 

(-^) f(^~) = (-j) I dx^ ■■■dx-e{x- -x:[)---e{x-_^-x-)f{x-) 
\0 / J-oo ^^_3Q^ 



1/1 1 
+ 



2\k+ + ie k+ -ie 
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Clearly the = modes are removed with this definition, the singularity of is regu- 
larized, and the vacuum is guaranteed to be trivial. Of course, the question now is what 
happens to the nontrivial QCD structure within a trivial vacuum? Unlike the sigma model, 
here it is not obvious that we can recover the effect of zero modes by the use of a shifting 
technique. However, although the fc+ = singularity is removed by Eq. (2.16), the infrared 
divergences from the small longitudinal momentum, i.e., surrounding the = region, 
are still present in the above Hamiltonian. These infrared divergences arc hidden in the 
three point quark-gluon and gluon-gluon vertices, and arc associated with the QCD vacuum 
properties. Recently, Wilson pointed out that, due to the different structure of the light 
front power counting, we can introduce noncanonical counterterms to remove the light-front 
infrared divergences [152] in light-front QCD, as we will discuss later. These counterterms 
may be the sources of the nontrivial vacuum effect. In order to understand the behavior of 
infrared divergences, we will first discuss how they occur in perturbation computations. 

V-3. Perturbative Light-Front QCD 

The light-front gauge singularity discussed above will lead to severe divergences even 
in perturbation theory, although Eq. (2.16) provides a well-defined regulator (a generalized 
principal value prescription, Eq. (5.30)) for the small k'^ momentum. In Feynman's per- 
turbation theory, the use of the principal value prescription still leads to "spurious" poles 
in the light-front Feynman integrals, which prohibit any continuation to Euclidean space 
(Wick rotation) and hence the use of standard power counting arguments for Feynman loop 
integrals [59]. This causes difficulties in addressing renormalization of QCD in Feynman 
perturbation theory with the light-front gauge. In the last decade there are many investiga- 
tions attempting to solve this problem. One excellent solution is given by Mandelstam and 
Leibbrandt, i.e., the Mandelstam-Leibbrandt (ML) prescription [57,58], which allows con- 
tinuation to Euclidean space and hence power counting. It has also been shown that, with 
the ML prescription, the multiplicative renormalization in the two-component light-front 
QCD Feynman formulation is restored [60]. 

Unfortunately, the ML prescription cannot be applied to equal-x^ quantization be- 
cause the ML prescription is defined by a boundary condition which depends on itself 
[151] and is not allowed in equal-x"'" canonical theory. Yet, as was pointed out recently 
by Wilson [152], light-front power counting differs completely from the power counting in 
equal-time quantization in that noncanonical counterterms are allowed in light-front field 
theory. In other words, multiplicative renormalization is not required in light-front QCD. 
Furthermore, the current attempts to understand nonperturbative QCD in light-front coor- 
dinates is based on the a;+-ordered (old-fashioned) diagrams in which no Feynman integral 
is involved. Thus the power counting criterion for Feynman loop integrals is no longer 
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available in light-front QCD Hamiltonian calculations. In Hamiltonian perturbation theory 
with the principle value prescription, light-front QCD contains severe linear and logarith- 
mic infrared divergences. Here I will give some results from the a;'*"-ordered perturbative 
loop calculations and renormalization of light-front QCD Hamiltonian theory up to one-loop 
[82,148], where the infrared divergences are systematically analyzed. Since light-front power 
counting allows noncanonical countcrtcrms, a complete understanding of renormalized light- 
front QCD may not be worked OT.it within perturbation theory; new renormalization and 
regularization approaches need to be developed, which wc will discTiss in the next section. 

The following calculations are performed in the formulation of the light-front Hamil- 
tonian perturbation theory with a two-component representation [82]. The light-front QCD 
Hamiltonian can be rewritten as a free term plus interactions, 

H = J dx-d^x±{no + Hint) = Hq + Hi. (5.31) 
It is easy to show that 

Wo = + (5.32) 

T^int = T^qqg + '^ggg + ^qqgg + ^qqqq + ^gggg ' (5.33) 



and 



-Hgggg = ^ f^' [aIAIA^^AI + 2^^) (A^5+A^) {Ap+ Ai) ^ 



(5.34) 



= 5</'^{-2 (^) {d-A^) + a- A^ (^^) {a -8^ + m) 
+ (-^) (a • 5j_ - m)a ■A_i_^(p, 

nggg = gr"'[d^A^M + {d^ K) (^) (^i^+^^j , (5-35) 

^qq99= 9^ ^V^^ O ■ A^(^-^^a ■ A^^ 

+2(^) (r^^A^a+A^j (^'^rv)} ^^-^^^ 

= '^qqggl + '^qqgg2 ) 

nqqqq = 2g' { ( ^) {^^ T'^ ^) ( ^) {^^ ^)\ , (5.37) 



(5.38) 



T^ggggl + 'Hgggg2 , 
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where is the hght-front fermion two- component representation (see Eq. (2.34) discussed 
in Sec. II). 

The x+-ordered Hght-front QCD perturbative theory can be obtained from the famil- 
iar perturbation expansion in quantum mechanics. The perturbation expansion of a bound 
state is given by (in the Rayleigh-Schrodinger perturbation theory) : 

I^) = E( f (5-39) 

where |$) is a unperturbative state, Q and are defined by: 

Q =!$)(*!, H'j = Hi-AE, AE = {^iHil"^) . (5.40) 

With this perturbative expansion formula, the mass, wave functions, and couphng constants 
renormahzations can be expressed as follows. For the convenience of practical calculations, 
we consider the expressions in momentum space. 

i). Wavefunction renormalization: In momentum space, the perturbative expansion 
of a state is given by 



I*) 



(l^) + E' 

I ni 



|ni)(ni|g^|$) 
P~ - Pni + ie 

(5.41) 



\m){n^\H'j\n2){n2\H',\<^ 



{P -Pni+ie){p -Pn2+i^) 



which has not been normalized, where |ni), |n2), • • • are properly symmetrized (antisym- 
metrized) states with respect to identical bosons (fermions) in the states and Yl' Eq. 
(5.41) sums over all intermediate states except the initial state |<I>). The normalized wave 
function is defined by 1^*') = ^/Z^\'^), where the factor Z$ is the wavefunction renormal- 
ization constant: 

= m = 1 + E't^^^^^^ + • • • • (5-42) 

m \P ~ Pni -r ie) 

ii). Mass renormalization. The mass correction can then be computed from the 
"energy-level" shift, i.e., the correction to the energy of an on-mass-shell particle. It is 
obvious that the perturbative correction to the light-front energy (p~) is given by 

Sp- = {^\{H - Ho)\^) = {^\Hi\^) 
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Using the mass-shell equation = p^p^ — p^, and recalling that p"*" and pj_ are the 
conserved light-front kinematical momenta, we obtain the mass renormalization in the old- 
fashioned perturbative light-front field theory: 

Sn? = p+5p- = p+imA^) + y+y' K^il-^_^l^)l' + . . . . (5.44) 

P -Pm+i€ 

iii). Coupling constant renormalization. The coupling constant renormalization is 
obtained by the perturbative calculation of various matrix elements of the vertices in Hj. 
Consider a vertex H} that is proportional to the coupling constant g, we have 

^ w {^f\H}\ni){m\H'j\^^) 
pT - Pni + «e 

^ y / {'i>f\H'j\ni){ni\H',\n2){n2\H}]\^i) (5.45) 
nt^2 (Pj ~ P'^^ + ^Pj ~ + ^e) 
^ y f {<i>f\H'j\ni){ni\H}\n2){n2\H'j]\^i) 

nt^2 (Pj ~ P^^ + ~ + 

^ ^/ (<^/|iJl|ni)(ni|g;|n2)(n2|g;]|a>.) ^ 

n'fea (^'?" - P'^^ + (^'i" - P^^ + 

where Zg is the multiplicative coupling constant renormalization, and and Zf are the 
wavefunction renormalization constants of the initial and final states. 

It is also convenient to express the above perturbation expansion in terms of the dia- 
grammatic approach. All matrix elements of H\ for the light-front QCD vertices arc listed 
in Table II with the corresponding diagrams. These are obtained by directly calculating 
the matrix elements between free particle states. The rules for writing the expression of 
perturbative expansions from diagrams are as follows: 

• Draw all topologically distinct x^-ordered diagrams. 

• For each internal line, sum over helicity and integrate using / for quarks 
and/^^^(A;+)forgluons. 



For each vertex, include a factor of 167r^(5^(p/ —pi) and a simple matrix element listed 
in Table I. 



• Include a factor (p^ — J2nPn + ^ l*^^ (Pf ~ J2nPn + ^1 fo^" ^ach intermediate 
state, where J2nPn sum over all on- mass-shell intermediate particle energies. 
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TABLE II. The a;"^-ordered Hamiltonian Diagrammatic Rules 
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• Add a symmetry factor S for each gluon loop coming from the symmetrized boson 



To illustrate the above computation scheme and to see the severe light-front infrared 
divergences, let me list some calculations up to one-loop based on the x^-ordered diagram- 
matical approach. 

i). Quark wavefunction and mass renormalization. Besides the infrared divergence, 
which is regularized by Eq. (5.30), there are also ultraviolet divergences for which we use 
a transverse cut-off regularization: |k±| < Aj^. For this simplest regularization scheme the 
one-loop light-front quark energy corrections (for the three diagrams in Fig. 5, respectively) 
are 



states. 



Spi = - 




2 2 

p — m 



2ln^--] InA^ 




(5.46) 
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FIG. 5. The a; "'"-ordered graphs for the one-loop correction to the quark mass and wave function 
renormahzation. 



where f{x) = {xm? —x{l — x)p'^). This shows that, in the one-loop quark energy correction, 
one-gluon exchange gives rise to both hnear and logarithmic infrared divergences. The 
instantaneous fermion interaction contribution (see in Fig. 5b) contains only one 
logarithmic divergence which cancels the logarithmic divergence in 5p{ . The instantaneous 
gluon interaction contribution (5p'^ of Fig. 5c) has a linear infrared divergence which 
precisely cancels the same divergence in Sp^ . This cancellation of linear infrared divergences 
is based on the use of the regularization for A:"*" ^ in Eq. (5.30) [82]. 

The quark mass correction (dropping the finite part) is then given by 

Sm^ = p'^Sp'\.„2^..„2 = -—^Cfrn^ln — ^, (5.48) 
^ 47r^ 

which is longitudinally infrared divergence free; and the quark wavefunction renormahzation 
constant is 

d5p~ 



Z2= 1 + 



dp 



The wavefunction renormahzation contains an additional type of divergence, the mixing of 
infrared and ultraviolet divergences, that does not occur in covariant calculations. This is 
the 'spurious' mixing associated with the gauge singularity. It corresponds to the so-called 
light-front double pole problem in the Feynman theory with the use of the light-front gauge 
and the principal value prescription that prohibits any continuation to Euclidean space 
and power counting in Feynman loop integrals. In the a; "'"-ordered Hamiltonian perturba- 
tive theory the power counting is different. The above argument of power counting for 
Feynman loop integrals may be irrelevant. Furthermore, since the second order correction 
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FIG. 6. The a;+-ordcred graphs for the one-loop correction to the gluon mass and wave function 
renormalization. 



to wavefunctions must be negative [see Eq. (5.42)], Eq. (5.49) shows that it is the additional 
infrared divergence that gives a consistent answer for wavefunction renormahzation. 

ii). Gluon wave function and mass correction. The one-loop light- front gluon energy 
corrections are given by (see Fig. 6) 



87r2 



■Ca 



21n- 



11 
~6" 



InA^ 



+ / dx(2- x(l -x)- -r^^ ln(x(l - x)q'^ 
Jo \ \x\ J 



+ 



Ai 



(5.50) 



^2 



flnAi 



— / dx{2x'^ — 2x + 1) ln(m^ — x(l — x)q' 
Jo 

21nAi-2 /^dxln(m2-x(l-x)g^) 
Jo 



m 



(5.51) 
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(5.52) 

oo 



where fc^ is the internal fermion momentum oo. In Eq. (5.50), there is also a 

mass singularity, which is well known in gauge theory and which is usually regularized by 
introducing a small gluon mass {q^ = iiq) in the energy denominators. The gluon mass 
and wavefunction renormalizations are 



2 ^^11 , (7+\„ . 2„ „ . A2 



l + { 21n-^ \Ca\^^ - —TfNflu^. 

87r2 \ V 6 e ) n% 3 ^ ^ m 



-2CAln'^| 



2 

(5.55) 



All severe divergences appear in the gluon sector: quadratic and logarithmic UV di- 
vergences, linear and logarithmic IR divergences, a gluon mass singularity and an unusual 
large longitudinal momentum logarithmic divergence. Only the linear infrared divergences 
are cancelled with the principal value prescription. The gluon mass correction is not zero. 
The first term in (5.54) is a fermion loop contribution (Fig. 6b), which is the same as 
the photon mass correction in QED. In addition, the gluon mass correction also contains a 
severe mixing of the quadratic UV divergences with the logarithmic IR and UV divergences 
of the longitudinal momentum. It is caused by the instantaneous fermion and gluon interac- 
tion contributions plus the tadpole effect of the normal four gluon interactions (Fig. 6c-6e) . 
This kind of divergence behaves in the same way from both the fermion contribution and 
the gluon contribution. The non-zero gluon mass correction of Eq. (5.54) is not surprising 
because it has the same divergence feature as the photon mass correction in light-front QED 
[Eq. (5.54) will be reduced to the photon mass correction when we set Tf = 1, Ca = 
and Nf = 1]. In a covariant calculation, the zero gluon mass correction is true only for 
dimensional regularization which "removes" or drops the mass correction. In the present 
calculation, maintaining zero gluon mass requires a mass counterterm, as is known in QED. 
The difference between QED and QCD is only manifest in the gauge boson wavefunction 
renormalization. For wavefunction renormalization, again there is an additional mixing of 
UV and IR divergences, which provides the correct sign for the wavefunction renormaliza- 
tion constant. Besides this feature, there is a contribution from the gluon loop (Fig. 6a). 
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As wc will see in the following, after the cancellation of the mixing divergences, it is this 
contribution that leads to asymptotic freedom in QCD. 

However, the above calculations show that the wave function renormalization contains 
several complicated pure infrared divergences and a mass singularity from the massless 
gluon. This complexity can be avoided if we introduce a mass scale fj, for the minimum 
cut-off for the transverse momentum: 

Ai > > (5.56) 

and assume that is much larger than all other masses in the theory. With this regulator, 
the quark and gluon mass and wavefunction renormalizations become simple, 

2 A 2 

Sm^ = -^Ctm^ In , (5.57) 



^{TfNfm^ In ^(Ai - /.^) - T,iV,) (l - In ^) } , (5.58) 



47r^ t II 



-2 ^ ^11 Q+\ 2^ 1 . A2 



1 + i hr - 21n ^ ) - ^TfNj (5.60) 



This shows that all the pure infrared divergences in the wavefunction renormalizations 
arc removed. The remaining unfamiliar divergences arc the mixing of quadratic UV with 
logarithmic IR divergences in the gluon mass correction, and the mixing of UV and IR 
logarithmic divergences in the wavefunction renormalizations. The mixing divergences in 
the gluon mass correction could be removed by a gluon mass counterterm, while the mixing 
divergences in the wavefunction renormalizations are cancelled in physical quantities. This 
can be seen from the coupling constant renormalization. 

iii). Coupling constant renormalization. For convenience, we set the external gluon 
momentum q{q^ , q'^J = 0. The quark-gluon vertex is then reduced: 

Vo = 2gT^,^6^,xA*- (5-61) 
In a; "'"-ordered perturbation theory, the one-loop vertex correction is given by 

SVo = {Vi + V2 + V3 + V4 + V5 + VejVo , (5.62) 
where V^, n = 1 — 6 are represented by Fig. 7. 
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FIG. 7. The a;+-ordered graphs for the one-loop correction to the quark-gluon vectex renormaliza- 
tion. 
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^3 = -i^ f 4 - 2 In ^) f-lc^ + Cf) In A , (5.65) 



^5 = , = . (5.67) 

To evaluate the contributions to the coupHng constant we have to multiply Vi and 
V2 by ^ in order to take into account the proper correction due to the renormalization of 
initial and final states [153,154]. Thus adding the contributions together, we have, 



SVo = (^^Vi + ^V2 + V2 + V4 + V5 + Ve]Vo 



(5.68) 



Note that all mixed divergences cancel now. The correction to the coupling constant is 
given by 

gR = g{l + 5g) = gi^l + -^[^^Ca- fiV/T/) In A | . (5.69) 

By redefining the bare coupling constant g such that g^ is finite. Thus we have given 
all renormalization quantities in QCD up to one-loop order based on the old-fashioned 
Hamiltonian perturbation theory. 

From these results, the anomalous dimensions for quarks and gluons and the (3 func- 
tion up to one-loop can be easily calculated. The anomalous dimension of the quark field 
to order g"^ is 

The momentum-dependent term implies that the quark anomalous dimension is gauge 
dependent. The anomalous dimension for the gluon field is 

which is also gauge-dependent. In the case of = 0, the gauge dependent term can 
be removed, and Eq. (5.71) is reduced to Gross and Wilczek's result in their Feynman 
calculation with = and g"*" = [42]. The (3 function is 



777 



778 



LIGHT-FRONT DYNAMICS AND LIGHT-FRONT QCD 



VOL. 32 



which is the well-known result to one loop order and is infrared divergence free, as we 
expected. 

Prom the above result, we see that there are severe light-front divergences in light- 
front QCD. Systematic control of these divergences is required a priori before we perform 
any practical numerical calculation in light-front coordinates for QCD bound states. Prom 

the basic one-loop calculations, one can see that, in the old-fashioned perturbation theory, 
light-front QCD involves various UV and IR divergences. Some of the divergences have not 
even been encountered in covariant and noncovariant Feynman calculations to the same 
order. Among various light-front divergences, there are two severe divergences one has 
to deal with in the old-fashioned theory for light-front QCD. The first is the mixing of 
UV and IR logarithmic divergences in wavefunction renormalization. The occurrence of 
the mixing divergences may not be a severe problem. The mixing divergences should be 
cancelled completely for physical quantities, as we have seen from the coupling constant 
renormalization. We expect that the problem of mixing divergences may not exist when we 
consider real physical processes. The second problem is the infinite gluon mass correction. 
In the old-fashioned Hamiltonian theory dimensional regularization is not available to avoid 
the nonzero gluon mass correction. To have a massless gluon in perturbation theory, we 
have to introduce a gluon mass counterterm. In the leading order (one-loop) calculation, 
there is no difficulty arising from a gluon mass counterterm. However, when we go to the 
next order, it has been found that the gluon mass counterterm leads to a noncancellation 
of infrared divergences. The non-vanishing infrared divergences could introduce non-local 
counterterms in both the longitudinal and transverse directions. In equal-time quantization, 
such non-local counterterms are forbidden for a renormalizable theory. Here, these non- 
local counterterms are allowed by the light-front power counting. This is a special feature 
of light-front QCD. One speculation from this property is that the non-local counterterms 
for infrared divergences may also provide a source for quark confinement [85]. 

In summary, renormalization in light-front QCD Hamiltonian theory is very different 
from conventional Feynman theory and it is an entirely new subject where investigations 
are still in their preliminary stage. In perturbative calculations, careful treatment could 
remove all severe infrared divergences for interesting physical quantities in light-front QCD. 
Por nonperturbative studies, the cancellation of severe infrared divergences may not work 
because certain approximations (e.g.. Pock space truncation) might be used. These ap- 
proximations may also break many important symmetries such as gauge invariance and 
rotational invariance. It is the hope of the current investigation of light-front renormaliza- 
tion theory that the counterterms for the light-front infrared divergences may restore the 
broken symmetries and also provide an effective confining light-front QCD Hamiltonian for 
hadronic bound states. 
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V-4. Infrared Divergences and Nontrivial QCD Properties 

Before we go on to discuss a possible approacli to tiie nonperturbative dynamics of 
light-front QCD, I would like to further address the infrared singularity and its relation 
with the nontrivial structure of QCD on the light-front [150]. 

The infrared divergences arise from the elimination of the unphysical gauge degrees 
of freedom. In coordinate space, infrared divergences are associated with the singularity at 
longitudinal infinity. For physical states, finite energy density requires that the longitudinal 
color electric field strength must vanish at the longitudinal boundary: 

^alx-=±oo = 0, (5.73) 

which is a condition to canonically remove the light-front infrared divergences and which 
leads to a constraint on the Aj^ at longitudinal infinity: 

= ^4 dx-j:{x-,x) . (5.74) 



2 

Eq. (5.74) is consistent with the choice of antisymmetric boundary condition from the 
definition of Eq. (2.16). Clearly, Eq. (5.74) is satisfied only for physical states. In per- 
turbation theory, we cannot use this condition because in perturbative QCD we consider 
not only physical states but also color non-singlet states for which Eq. (5.74) may not be 
satisfied. Therefore, the main effect of Eq. (5.74) should be manifested in nonperturbative 
dynamics, i.e., bound states. To address hadronic bound states, the existence of nontrivial 
structure, namely, confinement potentials and the mechanism of dynamical chiral symmetry 
breaking, is crucial. An explicit derivation of such nontrivial properties from QCD is still 
lacking. Here I will show how Eq. (5.74) can really infiuence the nontrivial properties of 
QCD. 

i). Topological winding number. Let us consider the axial current equation (for zero 
quark mass) 

« = ^/^Tr(F^.i?n, (5-75) 

where the axial current is = 'ipj^^5il^, and the dual field strength is F^^^ = ^e^'^'^T^p. 
The winding number in light-front QCD is defined as the net charge between x'^ = —oo 
and x"*" = oo, 

f 
87r2 

The integration on the r.h.s. of the above equation is defined in Minkowski space (M) and 
can be replaced by a surface integral. It has been found [150] that 

o _ 

AQ5 = -Nf-^ dx+d'^x±Tr{A-[A^,A^]) , (5.77) 

TT^ J x-=-oo 



^Q5 = Nf^l d^'xTriF^^rn- (5.76) 

M 
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where A'g^\x=±oo is determined by Eq. (5.27). We see that a non-vanishing AQ^ is gener- 
ated from the asymptotic fields of Al^,A~ and their antisymmetric boundary conditions at 
longitudinal infinity. 

For physical states, it is particularly interesting to see from Eq. (5.74) that the 
asymptotic physical gauge fields are generated by the color charge densities integrated over 
x". Thus, the topological winding number in the = gauge can be explicitly explored 
fromEq. (5.74). 

ii). Non-local interactions in the transverse direction. From Eq. (5.74) we see 
that the asymptotic A^^ fields at the longitudinal boundary are proportional to the color 
charge density in the transverse space and also that they involve non-local behavior in the 
transverse direction (induced by the transverse derivative). Intuitively, we may separate 
the transverse gauge potentials into a normal part plus a boundary part, 

Al = Al^ + Ais, (5.78) 

where 

^aN\x-=±oo = 0, d^^aN\x-=±oo = (5.79) 

5'4bU-=±oo = ^j{pl{x±) + Pl{x±)) . (5.80) 

In Eq. (5.80), Pa{x±) denotes the color charge densities integrated over x~ . The conditions 
ofEqs. (5.79) and (5.80) do not uniquely determine the separation of Eq. (5.78). Generally, 
there arc two types of separation for Eq. (5.78). One, is to consider A*^^ the long-distance 
fields generated by the boundary integrals and the short-distance fields determined by 
free theory. In this case, if we are only interested in the low-energy dynamics, the effect 
of the Al^j^ fields may be ignored. This separation is physically very interesting but it is 
practically very difficult to realise analytically. Another possibility is to choose a simple 
solution for the A^^^ that satisfies Eq. (5.80). In this case, the ^^tv have the trivial boundary 
condition, Eq. (5.79), but are not determined by the free theory. The Hamiltonian is then 
expressed only in terms of the ^^^v ^^'^ boundary behavior of transverse gauge fields 
is replaced by the effective interactions. A convenient choice for yl^^ which satisfies Eq. 
(5.80) is 

= -| / dx'-dx''e{x- - x'-)e{x' - x"){pi{x') + pl{x')) , i = l,2. (5.81) 

With such a formal solution, we obtain a new Hamiltonian in terms of A^^^ that contains 
many effective interactions induced by Eq. (5.74). All these effective interactions involve the 
color charge densities and involve non-local behavior in both the longitudinal and transverse 
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directions. One of the lowest order interactions, for example, is given by 



■oo 




(5.82) 



where r]^^ = 1(0) for i ^ (=)i- This is a linear interaction in both the longitudinal and 
the transverse directions. Hence, Eq. (5.74) leads to numerous many-body non-local color 
charge interactions which may lead to confinement. 

Still the Hamiltonian contains, in principle, an infinite number of many-body in- 
teractions generated by the boundary integrals (or obtained from the counterterms of the 
infrared divergences). This is a consequence of the boundary integrals in a non-abelian 
gauge theory due to the existence of nonlinear gluon interactions. It is also true in other 
gauge choices, such as the Coulomb gauge [155] or the axial gauge [157]. Practically, we 
do not know how to solve for hadrons from such a complicated QCD Hamiltonian. Never- 
theless, in the light-front QCD, one can set the vacuum to be trivial, the nontrivial QCD 
features for physical states are switched to the field operators and are manifested in the 
asymptotic behavior of physical gauge fields at longitudinal infinity. The trivial vacuum 
with nontrivial field variables in light-front QCD may provide a practical framework for 
describing hadrons. In the next two sections, we will explore these problems from a totally 
different point of view. 



As I mentioned in the very beginning, the theory of nonperturbative QCD (at the 
hadronic scale) has indeed not been defined at the present time. What I mean is that since 
we do not know how to renormalizc QCD nonperturbatively, all results we obtain (if we 
can calculate) may have nothing to do with the data observed in various experiments. In 
other words, the bare structure of QCD (or a field theory) does not directly connect to the 
real hadronic world although it is constructed from the beautiful requirement of various 
symmetries. Without knowing precisely the renormalization scheme for a field theory, the 
theory itself has not indeed been defined. 

One may argue that since we even do not know what is nonperturbative field the- 
ory, how can we talk about nonperturbative renormalization. It is true that there is no 
clear description in most field theory text books of nonperturbative field theory. However, 
quantum field theory is developed within the framework of quantum mechanics which is 
defined nonperturbatively in terms of the Hamiltonian formalism. From the descriptions 
in the previous sections, we have seen that light-front field theory can be formulated as an 
ordinary quantum mechanical system. This may provide a possible way to establish a non- 
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perturbative field theory for relativistic dynamics on the Hght-front, as we have emphasized 
throughout this paper. 

In this section, I will discuss a new renormalization scheme for Hamiltonian systems 
that was developed recently [83-85]. It begins with the light-front power counting plus a 
cutoff scheme on the constituents to construct the basic structure of the QCD Hamiltonian, 
then using a similarity renormalization approach to remove the cutoff dependence. The 
renormalizcd theory is given by an effective Hamiltonian that is characterized by a scale 
parameter. It demands that when the scale is taken at the hadronic energy level, the theory 
should produce all the physics of the hadrons. When the scale is far from the hadronic scale, 
the theory must reproduce all the results of perturbative QCD. 

VI- 1. Light-Front Power Counting: Canonical Structure 

Before discussing renormalization, I would like to introduce the light-front power 
counting rules which determine the possible structure of operators for the canonical Hamil- 
tonian [85]. 

Light-front power counting is done in terms of the longitudinal coordinate and the 
transverse coordinate x±. The power counting analysis is to deduce the most general struc- 
ture of divergences that arises from increasing the powers of the interaction Hamiltonian in 
perturbation theory. But power counting based on the kinematical symmetries of the light 
front is different from power counting based on the kinematical symmetries in equal-time 
coordinates. This is immediately transparent from the light-front dispersion relation for 
free particles, = -'-^^ — . Because the energy separates the fc+ and k± dependencies, the 
subtractions are not constants. For example, when k± gets very large the energy diverges 
no matter what A;+ is. Thus, in general, we get a divergent constant multiplied by a function 
of k'^ . In position space this translates into divergences at small x± being nonlocal in x~ 
and spread out over the light front. A similar result follows for the case when A;+ gets 
very small. It does scale differently under x~- and x_L-scaling (strictly speaking, a unique 
transverse scaling behavior holds only in the absence of masses). This situation is to be 
contrasted with the equal-time case where the relationship between energy and momentum 

in the equal-time theory is E = \/k^ + m? . In this equal-time form, if A;_l cx) while kz is 
fixed, the k^ dependence becomes negligible and arbitrary functions of kz cannot arise. 

i). Power counting analysis and canonical light- front QCD Hamiltonian. The scaling 
properties of the canonical Hamiltonian is determined as follows. Prom the canonical com- 
mutation relations (5.20) and (5.21), under the scale transformation 3^ ^ 3y Stic ij iZ/ _|_ ^ 

x'^ = txx, one can find the power assignments for the field variables 




(6.1) 
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The power assignments for the derivatives is obvious: 

, d+:^, -^-.x-. (6.2) 

From the Ught-front dispersion relation, we see that the Hamiltonian H = P~ scales just 
like x~ /x\, 

x" 1 1 
H : — > H : , m : , (6.3) 

CCj^ iCj^ X I 

where 7i and m are the Hamiltonian density and the mass parameter. Note that under 
scale transformations masses scale as constants. Thus the Hamiltonian does not have a 
unique scaling behavior when masses are present. Here the mass is assigned a power only 
based on dimensional analysis. 

The canonical light-front QCD Hamiltonian density is a polynomial in the six compo- 
nents m, A_i_, ^, 5*^, 9^, and The most general structure we can build for the canonical 
Hamiltonian density, which has dimension j^^, is 

(eCt)P(A^,ax,^)4-2p(5-h)-p. (6.4) 

Here the expression {Aj^,d±,m)^~'^^ stands for monomials of order 4— 2p in any combination 
of the three variables ^4^, d±,m, which immediately shows that the allowed values of p are 
0, 1, 2. The resulting canonical QCD Hamiltonian structure is 

p = 0: A\,Ald±,Aldl,'m?Al 

p = l: -^{(^^){Al,A±d±,dl,mA±,md±,m^), ^6.5) 

Physically, the terms that are proportional to the operator ^ correspond to the elimination 
of the dependent variables ip- and A~ . Specifying the precise way in which ^ acts, we 
can enumerate those terms which obey the canonical rules: 

p = 0: {A^)\ {A^fd^, {A^f{d^f,m^A^)^ , 

{A^d+A^)^{d^As.), {A^d+A^)(^^^\A^d+As_) , 

P= 1 : m[^)\A±d+A^),{CC^)^{d^A^), (6.6) 

im,d^,A^)e-^{{m,d^,A^)C}, 
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Comparing with the Hamiltonian derived in the last section, one sees that the free terms 
m?A\_ and m^'^ d_\_-^^ are absent. The absence of the first term is because of the presumed 
gauge invariance. A term like the latter does appear in the free part when '0_ is eliminated, 
but it is cancelled by a similar term and leads to chiral symmetry for the light-front free 
Hamiltonian. 

ii). Structure of light-front divergences. There are mainly two types of divergences 
in light-front field theory: ultraviolet and infrared divergences. The infrared divergences in 
light-front QCD arises from the elimination of the unphysical degrees of freedom. In the 
Hamiltonian formulation, as we have shown in the previous section, these divergences are 
manifested from the products of interaction Hamiltonians. Each term in the interaction 
Hamiltonian involves an integral over the product of three or four field operators, some 
of them carry high momentum and others low momentum. Power counting allows us to 
estimate the divergences arising from the high momentum operators after they are inte- 
grated out, and the remaining low momentum operators provide the operator structure of 
the corresponding countertcrms. 

Explicitly, let us consider some simple examples: Consider the second-order shift in 
the energy of a gluon coming from a two-gluon intermediate state. As candidates for Hj, 
we take 

Hii,) = ff/""^ J dz-dh^d'AiAlAi , (6.7) 

Hi(2) = -gr'" J dz-d^z^ ( aT^'^a) ^i^^^c • (6-8) 

The field operators are separated into low-momentum parts which contain wavelets of width 
5x±, 5x~ , and high-momentum parts, which contain wavelets of width Sy±, Sy~ . 

First let us consider the ultraviolet divergences. To avoid confusion with infrared 
divergences, we set 6x~ = Sy~ . For a candidate Hj we choose Hi(iy Remember that 
to produce an ultraviolet divergence at least two operators have to belong to the high- 
momentum sector, a). Consider ~ / dz~ (i^ Zj^{d'^ A^ A^)(y-^A-'^^^ and the energy shift 
AE iJjj-i^-gi^iJj^i), where E.D denotes the energy denominator. scales like j^^^i 
E.D produces a factor while the x wavelet operators in produce a factor 
(llj)^ and there is an overall factor from the integral {^j^)^ ■ Thus the scaling behavior 
of the energy shift is x {5y^)^ x x {^f x that is, AE « a 

quadratic ultraviolet divergence; b). Consider ~ / dz~ z^d"^ A^f^^^{A'' A^)(^y^ and the 

energy shift AE ~ H^i^^-j^Hj(^y The only difference from case (a) is that the x wavelet 
operators in produce a factor (Ifj)^" Hence the scaling behavior of the energy shift 

is jsir)^ ^ ^^y^f ^ TSyI7 ^ fe^^ ^ (fe^"^' ^^^^ ~ I&IF' ^ logarithmic ultraviolet 

divergence. 
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A complete analysis of light-front divergences and the corresponding counterterms 
based on light-front power counting and phase space cell approach is given in [85]. It has 
been found that the power counting rule for ultraviolet counterterms is the same (in 6x±) 
as for the canonical Hamiltonian that there can be at most four operators scaling as 
Thus the ultraviolet counterterms are built out of products up to fourth order in ^, , 
A± , and d± but they can have a completely arbitrary longitudinal structure, except for the 
restriction from the power counting and the longitudinal boost invariance. As a result of the 
constraint, the vertex counterterms can involve a priori unknown functions of longitudinal 
variables. 

Now consider infrared divergences. We set Sx± = dy±. Remember that to produce 
an infrared divergence at least two operators have to belong to the low sector. Thus, one 
may take ^ J dz-cPz^d'Al^{A'A^)^,), Hj^^) ~ J dz-d^z^{^d'A'^y^){A^d+A%y 

a). Consider the energy shift AE ^ Hj(^2)WD-^i('i)- ^H^) ^'^^^'^^ like 5y~ . The energy 
denominator produces a factor Another H^2) produces a factor 5y~ . Thus AE ~ 

dy^ , which results in a linear infrared divergence, b). Consider the energy shift AE ~ 
iJ/(2) Then, scales like 6y^ x (||^)- The energy denominator produces a 

factor Hj(2) produces a factor 6y~ . Thus AE Sx~ , which results in a logarithmic 
infrared divergence. 

In general, the counterterms for infrared divergences involve arbitrary numbers of 
quark and gluon operators. They have a complex nonlocality in the transverse variables. 
This is in contrast to the divergent counterterm for the canonical instantaneous four-fermion 
interaction which is local in the transverse direction. If we identify the counterterms arising 
from infrared gluons (small longitudinal momentum) as the source of transverse confine- 
ment, then the unknown nonlocal transverse behavior would have to include confining effects 
at large transverse separation. 

One might worry that the appearance of functions of momenta in the counterterms 
could destroy the predictive power of the theory and lead to non-renormalizability. However, 
these functions are needed to restore Lorentz covariance and gauge invariance as g ^ Qs 
in a weak coupling treatment (see the discussion in the next section), and without such 
counterterms physical quantities will not even approach finite limits as cutoffs are removed 
for any coupling. 

VI-2. A New Cutoff Scheme for the Light-Front Hamiltonian 

The light-front infrared and ultraviolet divergences discussed above should be regu- 
larized properly in order to develop a renormalization scheme. To provide a renormalization 
process for a realization of nonperturbative Hamiltonian computation, we need to develop a 
cutoff procedure that is applicable to the Hamiltonian as a whole rather than just to order 
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by order perturbativc calculations that we have used in the previous section. One way to 
accomplish this is to cut off the single particle momenta appearing in the field variables 
themselves. Cutoffs on the constituent momenta A;+ and k± are called "constituent cutoffs" 
[85]. The resulting cutoff boundary on each constituent's relative momenta {k^,ki±) is 
given by 

kl = 2A^4F(^-4F)-m' for ^Po^ < < ^ , (6.9) 

and 



kL = ^-m^ for ^<kt, (6-10) 

where m is the lowest of the constituent masses. In order to define such a constituent 
cutoff which still limits the invariant mass of a state, however, we are forced to introduce 
a longitudinal momentum scale P^. Dependence on this scale must also be eliminated as 
part of the renormalization process. Then, in a considerable range of center of mass domain 
(Pj_ = 0, P"*" = to Pj_ = P+ > "x"), when the above cutoff constituent appears 
in a state, the internal mass of the state is guaranteed to be at least A^. 

For future practical numerical computations, we will also provide a buffer zone outside 
this constituent momentum cutoff boundary extending kf^ roughly by a factor of 2 and k^ 
by a factor of 1/2. We can accomplish this by setting the outside of the buffer zone at 

kL = 4A^^ (l - ^) - <kt<^, (6.11) 

and 



p+ 

kf^ =k^ -rr? for Kkf . (6.12) 

The buffer zone allows the use of a smooth cutoff in order to let the interactions die grad- 
ually. The cutoff A eliminates both ultraviolet transverse and infrared longitudinal degrees 
of freedom. 

Practically, the cutoff on constituents means that these cutoff boundaries are directly 
employed in the integrals over momenta in the momentum-space expansion of the field 
operators. Thus, the quark field operator is 

^{x) = E j[ '^^l^l^^ Mk)e-^'-^ + 4 A(fc)e-^-] , (6.13) 
and the gluon field operator is 

^'(^) = E X 167r3fc+ + 4*«i(^)e^''1 , (6-14) 
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where color indices have been suppressed. If we choose sharp momentum cutoffs, the 
momentum integrals may be explicitly written as 



dk+cPk± = jj dk+ j d^k±e{2A'^x{l - x) - ni^ - kl) 



(6-15) 

where now k^^^ = j^Po~- It should be noted that sometimes a sharp cutoff may cause 
a so-called non-analytic divergence in a practical calculation [85]. To avoid such non- 
analyticities in the structure of counterterms, a smooth cutoff can be introduced on the 
momentum integrals so that 

where x = k'^ /Pq . C{y) is a cutoff function which equals 1 for y = and decreases 
analytically to for y = 1. The integration over longitudinal momentum is thus cutoff at 
a minimum fc^^j^^ = ■ 

VI-3. Sigma Model with and Without Zero Modes 

Using the power counting argument and a cutoff scheme, we may generate a theory 
(even at tree level) that may have already manifested some nontrivial properties with a 
trivial vacuum. Here we consider again the sigma model to show such a possibility. Ob- 
viously, with a cutoff on the constituents, the zero modes are removed and the vacuum is 
trivial. To construct the effective Hamiltonian we can rely on power counting and locality. 
The Hamiltonian can generally be written as 



j dx <^x± 



-^(5x0 • -I- (?j_7r • 5j_7r) -I- V 



(6.17) 



Since the symmetry is broken only in the cr sector, V should be even in the 7r-field so that 
the symmetry under tt — ^ — tt remains. Also, since no inverse powers of mass are allowed, 
we can write 

y = y^'/'^ + + Ak/-^ + A27r^ + Ag^V^ + U<p^ + Asvr^,/. . (6.18) 

On the other hand, the current operator can also be constructed by power counting. 
The charge Q (= ^ / dx^ (f'x±J~^) is dimensionless, which leads to the canonical dimension 
of J"*" being ~ '^J^^y- Meanwhile, each term in d^J^ (= \{d^ + d^J^) — ■ J^) 
should have the same dimensions, which determines the dimensions of the other components: 
~ /..'*' ) J ~ r^TT- The components of J'^ are to be constructed from the operators 
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TT, (/>, S"*", d~ , and d'^ and the constants. The constants arc allowed to have the dmiensions 
of a negative power of x± (as masses do) but no power of . Since the operator ^ does 
not appear in the canonical scalar field theory, the operator ^ is not allowed to appear in 
the canonical current operator. One should also implement the symmetry that has to 
change sign under tt —it. Then the allowed structure of the components of is 

J+ = aiS'^TT + a2(9+7r)^ + a3(5"'"0)7r , 

= hldx.TT + h2{d^(P)TT + b^{dx.TT)(t), (6.19) 
J~ = Cld~ TT + C2{d^ TT)(t) + C^{d^ (t))TT . 

Now all of the coefficients in the current and the Hamiltonian can be determined by the 
requirement of current conservation, d^j^ = 0, and covariance. The results are: 

oi = 6i = ci, a2 = -62 = C2 = -as = 63 = -C3 = 1 , 

1 A (^-20) 

Ai = A2 = yAs = — , A5 = Aai, A4 = Aai, m\ = ml + 2\a\ , 

and 

ai = or m^r = . (6.21) 

If we take oi = then the potential is reduced to the canonical form, and the current is 
also of the canonical form. This corresponds to the full canonical theory with a symmetry 
preserving vacuum and a doublet in the spectrum. If we choose rriT^ = 0, we have 



(6.22) 

V 

which is the theory with zero modes eliminated and the pion is masslcss; then Hamiltonian 
explicitly breaks the symmetry, the vacuum is trivial, and there is no longer the notion 
of spontaneous symmetry breaking. Current conservation is preserved and vastly reduces 
the number of free parameters present in the effective Hamiltonian constructed from power 
counting and forces the pion to remain massless. The results are the same as we have 
obtained by the use of the shift approach in Sec. IV. We have argTicd that the shift 
technique cannot be applied to QCD. But the above solution is purely based on light-front 
power counting and some symmetry consideration, and it is generic for various theories. Of 
course, it should be emphasized that the sigma model discussed above is only at the tree 
level. For QCD, one must consider radiative corrections and the subsequent renormalization 
effects, which are much much more difficult. 
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VI-4. Similarity Renormalization Scheme 

Now, I am going to introduce a new renormalization procedure for a Hamiltonian 
formulation that was recently developed by Glazek and Wilson [83,85], which will allow us 
to construct an effective Hamiltonian for light-front QCD. The previous discussion shows 
that we can write down a light-front QCD Hamiltonian based on the power counting ar- 
gument and some requirement of symmetries, with a cutoff on constituents. To begin with 
such a bare cutoff Hamiltonian and to approach the hadronic solution we need to develop 
a renormalization procedure. The rcnormaUzation should allow us to remove all cutoff de- 
pendence, yet this is only a minor requirement. The highlight of the renormalization we 
must develop is that it allows us to reduce the effects of relativistic momenta in an effec- 
tive Hamiltonian which is dominated in weak coupling by nonrelativistic relative momenta. 
Then, the bound states of the effective Hamiltonian in weak coupling are pure qq, pure qqq, 
or pure gg states, exactly as predicted in the CQM. All major effects of gluon emission and 
absorption — or more complex processes — are absorbed into the effective Hamiltonian for 
these states. Once such an effective Hamiltonian is constructed, the computation of bound 
states involves readily executed numerical computations. 

The renormalization procedure discussed next to realize these requirements is called 
the similarity renormalization scheme for a Hamiltonian formulation. The basic idea is to 
develop a sequence of infinitesimal unitary transformations that transform the initial bare 
Hamiltonian Hb to an effective Hamiltonian H^^, where a is an arbitrarily chosen energy 
scale. 



The basic goal for the transformation is that H„ should be in band-diagonal form relative 
to the scale a, namely, the matrix elements of involving energy jumps much larger than 
a (other than jumps between two large but nearby energies) will all be zero, while matrix 
elements involving smaller jumps or two nearby energies remain in H^. The similarity 
transformation should satisfy the condition that for a — >^ oo, Her Hb and ^ l- The 
effective Hamiltonian we seek involves H^r with a on the order of the quark or gluon mass. 
Consider an infinitesimal transformation, then Eq. (6.23) is reduced to 



which is subject to the boundary condition limCT^oo ^ct = Hb- To meet our demand, we 
need to specify the action of T^. This can be done by introducing a scale a: 



Hfj = S^HbSI- . 



(6.23) 



dH^ 



(6.24) 



da 



X, 



Ei — Ej 



(6.25) 



(713 ~ 



Ei + Ej + a 
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into a function /o-jj = f{xaij) such that 

< X < -J, f{x) = 1 {near diagonal region); 

-J < X < /(x)drops from 1 to (transition region); (6.26) 

-|-<x<l, f{x) = (far off diagonal region); 

where f{x) is to be infinitely differentiable throughout < x < 1, including the transition 
points 1/3 and 2/3. Then we rewrite Eq. (6.23) as follows: 

dH„ii „ „ , Mr, f \H 

' {^'-'^ J crij ) ^ crij ) 

(6.27) 



— f aij\H ,T(j\ij -\- (In ff^ij) H(^^j , 



'^crij — ^. fijij)[Hla,Tu]ij (In /g-jj )i?cr«j • 



Since /(x) vanishes when |x| > 2/3, one can see that Haij does indeed vanish in the 
far off-diagonal region which will help identify divergent terms and determine the form of 
counterterms necessary to remove these divergences. It also can be seen that T^ij is zero 
in the near-diagonal region which means that a perturbative solution to Huij in terms of 
Hf-j will never involve vanishing energy denominator s.^ 

The solution for Hj^^ and Tg- are 

Hi^ = Hi + [Hi,. , T,.]^ , = H^^ + [Hi.',T,'] ^ , (6.28) 



where Hf^- = fmjHf- and the linear operation R is 



/•oo 

-faij / da'X^'ij , (6.29) 



Finally, one can obtain an iterated solution for Hi^j where all higher order terms have the 
substitution • • • , 

= Hl+[Hl,,Hl,rp] 

" 'R 

- . (6.31) 

+ [Hfa' , [Hf^t , Hlnj] ] H 

" ^ 'R 

® As we know, a vanishing energy denominator in the old-fashioned Hamiltonian perturbative calculations 
is a potential source of large errors in other perturbative renormalization schemes. 
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The counterterms in Hf must be chosen to cancel the divergences which occur in integrals 
over intermediate states at higher orders in Hf, and such counterterms must also then 
be included in higher-order iterations. If a limit to this process exists the Hamiltonian is 
said to be renormalizable. Now we see that the effective hamiltonian H^r causes transitions 
only between states staying "close to the diagonal" due to the factor f^ij, and the effects 
of transitions to and from intermediate states which are "far off-diagonal" have explicitly 
appeared in the effective Hamiltonian H^j as a perturbative expansion. 

Thus, using the similarity renornialization scheme, one transforms the Hamiltonian 
Hb into a manageable, band-diagonal form H^j. Hb is the bare cutoff Hamiltonian, forced 
to be finite by the imposition of some cutoff A. After we determine the form of the countert- 
erms and include them in Hb, each of the matrix elements of the transformed Hamiltonian 
Huij has no large dependence on A. Thus as we send A — > oo, we get a renormalized, 
scale-dependent effective Hamiltonian H^. This does not finish the renormalization of the 
Hamiltonian, however, for the finite parts of the counterterms in Hb will produce in H^ 
unknown constants and functions of momenta which must be adjusted to reproduce phys- 
ical observables and to restore the symmetries which were broken by the cutoff A. These 
quantities are to be fixed by solving H^; one should be able to do this with a combination 
of few-body Hamiltonian methods and weak-coupling diagrams. 

One can use the similarity transformation to bring the effective Hamiltonian to any 
scale (T„ and obtain a set of sequences {i?^}. As we change a we change the characteristic 
scale, but the physics is invariant with respect to this change for large enough A^v- The 
utility of the similarity renormalization scheme is that the transformation is invertible. 
Thus, if one finds a Hamiltonian that is finite and A- independent for any one scale a"„, the 
differential similarity framework guarantees that one can obtain a Hamiltonian that is finite 
and cutoff independent for all a. Note that we require that each matrix element of H„ be 
cutoff independent for external momenta which are small in comparison to the cutoff, which 
is more restrictive than just requiring this of the eigenvalues of iJ^- 

VII. A WEAK COUPLING TREATMENT OF NONPERTURBATIVE QCD 

In this section, I will discuss a new approach to QCD that was recently proposed based 
on the theoretically successful approach to bound states of QED and the phenomenologically 
successful approach of the CQM [85] . The basic question we try to answer is whether we can 
set up QCD to be renormalized and solved by the same techniques that solve QED: namely, 
a combination of weak-coupling perturbation theory and many-body quantum mechanics, 
from which we can, at least, reproduce the phenomenological success of CQM for hadrons. 
The renormalization scheme discussed in the last section may provide us with a practical 
framework to apply these two approaches together to QCD on the light-front. 
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To develop such a framework, one needs to change from the standard approach in 
QED by: a) the use of hght-front dynamics; b) the use of nonzero masses for both quarks and 
gluons; c) the need to handle relativistic effects which give rise to, for example, asymptotic 
freedom in QCD, which in turn leads to a fairly strong renormalized coupling constant and 
hence relativistic binding energies; d) the presence of artificial stabilizing and confining 
potentials which vanish at relativistic values of the coupling constant but nowhere else; 
and e) the suitable treatment of light-front longitudinal infrared divergences, which cancel 
pcrturbatively in light-front QED because of gauge invariance, and which is the path for 
seeking the source of artificial potentials in light-front QCD. 

The formulation starts with nonzero masses for both quarks and gluons from which 
one can consider an arbitrary coupling constant g which is small even at the quark-gluon 
mass scale. In the beginning, this will sacrifice manifest gauge invariance and Lorentz 
covariance. These symmetries are only implicitly restored (if at all) when the renormalized 
coupling is increased to its relativistic value, which we call Qs- The value Qs is a fixed 
number measured at the hadron mass scale. For smaller values of g the formulation lacks 
full covariance and is not expected to have the predictive power of QCD, but it allows 
phenomenology to guide renormalization and is defined to maximize the ease of perturbative 
computations and extrapolation to gg- The nontrivial QCD structures arise from non- 
cancelling divergences in the new framework, which are necessarily the sources of true 
confining potentials and chiral symmetry breaking in QCD. With the cutoff Hamiltonian 
one can have a trivial vacuum. We expect light-front infrared divergences to be sources of 
confinement and chiral symmetry breaking because these arc both vacuum effects in QCD 
and vacuum effects can enter the light-front theory through light-front longitudinal infrared 
effects, as we have explored in the previous sections. Because of the unconventional scaling 
properties of the light-front Hamiltonian, these effects include renormalization counterterms 
with whole functions to be determined by the renormalization process. 

The basic motivation of the approach is physical rather than mathematical. Phys- 
ically, a Hamiltonian with nonzero (constituent) quark and gluon masses and confining 
potentials is closer to the physics of strong interactions than a Hamiltonian with zero mass 
constituents and no confining potentials. Then, in the spirit of QED, renormalization ef- 
fects analyzed with the confining potential itself can be treated perturbatively, but only for 
generating an effective few-body Hamiltonian which can be solved nonperturbatively. In 
this section, I will discuss the scheme of how to address these arguments. For a detailed 
derivation, please look at our recently published paper [85]. 

VII- 1. Bare Cutoff Hamiltonian with Artificial Interactions 

We begin with a bare cutoff light-front QCD Hamiltonian in which the quark and 
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gluon are massive constituents, 

H'^iA) = J dx-dW^f^d^Aid^Ai + mlAf + (^^^^) ^] + Hf, , (7.1) 

where Hf^ is the interaction part of the canonical light-front QCD Hamiltonian given in 
Sec. V, the only difference with the superscript B is that all constituent field operators are 
imposed by a cutoff of Eq. (6.13) or (6.14) on their momentum expansion. 

The choice of light-front dynamics, massive quarks and gluons, and a particular 
cutoff scheme eliminates the traditional barriers to a weak-coupling treatment of QCD and 
allows us to begin with a trivial vacuum. In the equal-time theory, the QCD vacuum is 
thought to be a complicated medium which presumably provides both confinement and the 
spontaneous breaking of chiral symmetry. Here the vacuum is trivial so we have to find other 
sources for confinement and spontaneous chiral symmetry breaking in the cutoff theory. In 
the bare cutoff canonical Hamiltonian, one particular term of interest is the instantaneous 
interaction in the longitudinal direction between color charge densities, which provides a 
potential which is linear in the longitudinal separation between two constituents that have 
the same transverse positions (see the canonical light-front QCD Hamiltonian in Sec. V). 
In the absence of a gluon mass term, this interaction is precisely cancelled by the emission 
and absorption of longitudinal gluons [82] . 

The presence of a gluon mass term automatically prevents unbounded growth of 
the running coupling constant below the gluon mass scale and provides kinematic barriers 
to unlimited gluon emission. It eliminates any equal-time type infrared problems. The 
cutoff procedure on constituents for the Hamiltonian developed before is valid also only 
with nonzero quark and gluon masses in a specific frame, while a large number of states 
(the upper limit of their invariant masses is guaranteed to be above a large cutoff) are still 
available for study even in the boosted frames (see the discussion on the cutoff scheme in 
the last section). However, the nonzero gluon mass results in the cancellation of the linear 
longitudinal interaction being incomplete. The canonical Hamiltonian combined with a 
one-gluon exchange term (in momentum space) is (see Fig. 8): 



FIG. 8. The a;+-ordered qq interaction via one-gluon exchange. 
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from which one can find an additional hnear interaction due to the nonzero gluon mass, 
which is proportional to 



g2rparpa 



1 



ml 



(7.4) 



(Pl - P3 )^ + (Pl-L - P2±)^ ' 

Yet, this interaction falls off too rapidly in the transverse direction, but is too strong in the 
longitudinal direction, which leads to an instability even at weak-coupling. To remove the 
above instabilities in hadronic bound states, we need to subtract this linear interaction, 

— / (fx±dx dy jt{x ,x±)\x - y~\jt{y~ ,x±) , (7.5) 

from the bare cutoff hamiltonian. 

The cutoffs will violate Lorentz and gauge symmetries, forcing the bare Hamiltonian 
to contain a larger than normal set of counterterms to ensure a finite limit as the cutoffs 
are removed. Of special interest are counterterms that reflect consequences of zero modes 
(namely, modes with k'^ = 0), and thereby the effect of confinement and chiral symmetry 
breaking in the full theory. According to power-counting arguments, the counterterms for 
longitudinal light-front infrared divergences may contain functions of transverse momenta; 
there exists the possibility that the a priori unknown functions in the finite parts of these 
counterterms will include confining interactions in the transverse direction. The g ^ Qs 
limit may then be smooth if such confining functions are actually required to restore full 
covariance to the theory. 
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The task of solving the hght-front Hamiltonian at the relativistic value of g is far too 
difficult to attempt at the present time. Instead, it may be helpful to define a plausible 
sequence of simpler computations that can build a knowledge base which enables studies 
of the full light-front Hamiltonian to be fruitful at some future date. A crucial step in 
simplifying the computation is the introduction of artificial potentials in the Hamiltonian. 

The basic need is to incorporate the qualitative phenomenology of QCD bound states 
into the artificial potential. This qualitative phenomenology comes from three sources: ki- 
netic energy, Coulomb-likc potentials, and linear potentials. Three terms should be present 
in the weak-coupling Hamiltonian and all should have the same overall scaling behavior 
with g in bound state computations, namely g^, just as QED bound state energies scale as 
e^. The Coulomb- like terms can be constructed directly from Eq. (7.2). By taking a non- 
relativistic limit, the above effective interaction in the qq sector due to massless one-gluon 
exchange is reduced to 

2 1 
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By using an interpolating Fourier transformation. 
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can further be expressed in position space with longitudinal boost invariance: 



Vc = 9' 



+ 



m2c 



Pi ri Pi r2 



This is the light-front Coulomb potential, where r is defined as 



r = 



(7.8) 



(7.9) 



Sx~ is the light-front longitudinal separation of two constituents, 6x± the transverse sep- 
aration, rric the constituent mass, and p'^ the constituent longitudinal momentum. The 
p'^/rric in the definition of r ensures that the dimensions match. Of course, the positive 
or negative SU (3) charges must also be inserted. It is easy to check that the light-front 
radial coordinate is invariant under longitudinal boosts. However, it is not invariant under 
transverse boosts. Relativistically, the spinor matrix elements are different in different he- 
licity sectors. In the nonrelativistic limit this helicity dependence vanishes and we get the 
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same interaction in all helicity sectors. In practical calculation, these helicity dependent 
interactions in Eq. (7.2) must also be included. 

Finally, we need a linear potential term that is proportional to r and r' . In Coulomb 
bound states both are of order l/g^. Hence, to achieve an energy of order g^^ the linear 
potential must have a coefficient of order g^. Thus the linear potential term, in position 
space, would be proportional to g^r. To be precise and get the dimensions straight, its form 
is 



Vl = i(i 



rrvir m'r 



p+ p'+ 



(7.10) 



with P a numerical constant. 

The linear potential needs to exist between all possible pairs of constituents: qq, qq, 
qq, qg, qg, and gg, where q, q and g stand for quark, antiquark and gluon respectively. The 
linear potential must always be positive (confining) rather than negative (destabilizing). It 
cannot involve products of SU (3) charges as the Coulomb term does. The Coulomb term 
could be given a Yukawa structure rather than the pure g'^ /r term due to the nonzero gluon 
mass. All of the potential have to be Fourier transformed to momentum space and then 
be restricted to the allowed range of both longitudinal and transverse momenta after all 
cutoffs have been imposed. 

A primary rule for the artificial potential is that it should vanish in the relativistic 
limit. For example, the artificial potential might have an overall factor (1 — g"^ /g1) to 
ensure its vanishing at g = gg. This rule leaves total flexibility in the choice of the artificial 
potential since no relativistic physics is affected by it. To ensure that the artificial potential 
vanishes at gs without destroying its weak-coupling features, the subtraction term, Eq. 
(7.5), of the longitudinal linear potential has to be treated with care. It is suggested that 
the subtracted linear potential be multiplied by (1— ^^/^f) so that the subtraction begins to 
be negated only in order g^ , which is smaller for small g than the artificial g^ linear potential 
that needs to be dominant. All other terms in the artificial potential can be multiplied by 
(1 — g'^ / g1) instead. To ensure that the Coulomb term shows Coulomb behavior, at least 
roughly, at typical bound state sizes, it is important that the mass used in any Yukawa-type 
modification of the Coulomb term scale as g'^ rather than being a constant mass. 

The artificial potential is to give the weak-coupling theory a structure close to the 
CQM so that past experience with the quark model can be used to determine the precise 
form of this potential and to fit it to experimental data. Thus, an initial calculation 
involves QCD complications only in a very minimal form. The artificial potential is also to 
incorporate a linear potential in both the longitudinal and transverse directions to ensure 
quark confinement for any g. This is important for phenomenology. 

It is necessary to study where the artificial linear potential might originate from 
the counterterms of infrared divergences, and, especially in the relativistic limit, where 
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the artificial potential vanishes. The light-front infrared singularities give rise to both 
linear and logarithmic divergences. The linear divergences, however, contain the inverse 
of the longitudinal cutoff ^, which violates longitudinal boost invariance; hence, only the 
logarithmic infrared divergences can be the source of the linear confinement interaction. We 
expect that the logarithmic infrared divergences could occur in the 5^-order correction to 
Fig. 8, from which the counterterms might determine the structure of such linear potentials. 
But so far this calculation has not been completed. 

Now, we add the artificial potential to the cutoff canonical Hamiltonian. 



HeiA) = J dx-(fx^^l^d'Aid'Ai + mlAf + 



-V^ + ml 



d+ ri (7.11) 



where 



■\^-^j+{x)Vc{x,y)j+{y) + j+(x)VL(x,y)j+(y)| (7.12) 

'^{^~^^ J dx~d^x±dy~j+{x~,x±)\x^ -y~\j+{y~,x±) , 

j+ is the color vector charge density and the color singlet charge density. This start- 
ing bare cutoff Hamiltonian will allow us to determines an effective Hamiltonian via the 
similarity renormalization scheme. 

VII-2. Effective Light-Front QCD Hamiltonian and Bound State computations 

In this subsection we outline the construction of the effective light-front QCD Hamilto- 
nian for low-energy hadrons, and discuss the increasing levels of complexity in the light-front 
QCD bound state computations in this formulation. 

In the similarity renormalization scheme, the effective light-front QCD Hamiltonian 
that will be used to compute the hadronic bound states is 



— lim faij 

A—^oo 



ttB I ttB ttB ( 9<Jijk dajik \ 



(7.13) 



k ^ j k ^ i ^ k 

Here the bare cutoff Hamiltonian is divided as Hb = Hq + Hf . For the determination of 
the effective Hamiltonian , the unperturbed part of the bare Hamiltonian Hq is chosen to 
be that of free massive quarks and gluons with the standard relativistic dispersion relation 
in light-front kinematics. The interaction part Hf then contains the canonical interaction 
terms Hf, the artificial potential plus a complete set of counterterms. Recall that Hf^ 
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has the cutoff A which violates both longitudinal and transverse boost invariance. First 
one needs to identify the counterterms that must be included in so that the matrix 
elements of the effective Hamiltonian H^ij have no divergent dependence on A. Then one 
can send A — >^ oo so that the cutoff dependence is removed from the effective Hamiltonian 
Hfj. As in standard perturbative renormalization theory, A dependence will be removed 
order by order in g^, where g^^ is the running coupling constant at the similarity scale in 
the Hamiltonian matrix elements. This avoids having to solve the nonperturbative bound 
state problem in order to identify and remove the A dependence. 

Consider the effective Hamiltonian H(j generated to second order in the renormalized 
coupling constant. The counterterms H^^ up to this order contain i). the counterterms 
for the canonical instantaneous interactions, which are 



and ii) . the counterterms from one-gluon exchange given in momentum space, plus a similar 
term in the gluon-gluon sector. There are also quark and gluon mass counterterms, some of 
which have given explicitly in [85] . The effect of these counterterms is to completely cancel 
the leading radiative corrections from instantaneous gluon exchange. 

There is a question whether H„ij satisfies boost invariance after A ^ oo. There can 
be finite terms in which violate boost invariance yet cannot justifiably be subtracted. 
In this case violations of boost invariance can only disappear at special values of g where 
the coefficient of the boost violating terms vanishes. Clearly, one such special value has to 
he g = gs, as part of the restoration of covariance at gs- In the following it is assumed that 
the effective Hamiltonian has boost invariance simply as a result of taking the oo 
limit. is generated as the first step in solving the bare cutoff Hamiltonian Hb- 

The second step is to construct bound states from which are determined by the 
bound state equation of Eq. (4.11) with i^^- as the Hamiltonian. One can solve this field 
theoretic bound state problem in the standard fashion, using the bound state approach of 
QED. This requires us to identify a part of i/^j -ffaO) which is treated nonperturbatively to 
produce bound states. The essential simplification that makes further calculation possible 
is that H^Q does not contain any interactions that change particle number, so that the 
methods of few-body quantum mechanics can be used to solve this initial nonperturbative 
problem. All field theoretic corrections the interaction part of the effective Hamiltonian 
that arise from particle creation and annihilation will then be treated perturbatively, as in 



Now, one might ask, what is the dependence of this effective Hamiltonian H^j on 
the scale a and what is the physical consequence of the scale dependence? Suppose we 
restrict the states to only a qq pair. Then only number conserving (potential-like) terms 




(7.14) 



QED. 
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in H(j contribute. Now, if wc let the scale a approach a very large number, only the orig- 
inal interactions in Hb survive since the factor (1 — faikfajk) 0- Thus the canonical 
four-fermion interaction together with all counterterms survive, whereas the transversely 
smeared four-fermion interaction diminishes in strength. In the limit cr oo we are recov- 
ering Hb, which is no surprise. In this limit the effects of higher Fock sectors (for example 
qqg) can be recovered only by including them explicitly, and so it clearly becomes a poor 
approximation to include only the qq sector when a becomes too large. 

By lowering the similarity scale a, we reduce the allowed range of gluon momenta 
which can contribute, for example, to the binding of a quark and antiquark in a meson. 
These effects must appear elsewhere; we see that, through the similarity transformation, 
they are added directly to the Hamiltonian via the second and higher order terms in Hb- 
Thus, by lowering a, we put the bare gluon exchange effects of Hb into a qq potential in the 
effective Hamiltonian H^ perturbatively. This clearly changes the character of the bound 
state calculation. It changes from a field theoretic computation with arbitrary numbers 
of constituents to a computation dominated by an effective qq potential. If we choose the 
similarity scale a to be just above the hadronic mass scale the major effects come from 
the qq sector. The resulting nonrelativistic calculation will not see the scale a in the first 
approximation, since only states close to the diagonal, for which ff^ij = 1, will contribute. 
This can be outlined with the following schematic equation: 

H^CD = H^{9s)+Hf{g,) 

= {H^{gs) + H^,uf^c^al{9s)] + {Hf{gs) - H^rtificialigs)] (7.15) 

= H^o{g) + H^i{g) . 

Equivalents, H^o{g) = H^{a,g^) + (1 - gl/g'^,)H^^^ifi^i^i{a, g^) and H^i{g) = Hf{a,ga). 

A basic goal of this approach is to construct a sequence of computations in light- 
front QCD with growing levels of complexity. The major sources of the complexity are 
analyzed in [85]. The further computations are i). an initial nonperturbative calculation 
of the hadronic masses to determine the parameters in H^q in the hadronic scale with the 
binding energies up to the fourth order of the coupling constant. This is a straightforward 
calculation; ii). the next, the leading radiative correction to these masses, which requires a 
calculation of the effective Hamiltonian up to g'^ and the binding energies up to g^, in which 
logarithmic infrared divergences should occur and the transverse linear potentials as the 
counterterms of these divergences could be justified; iii) the extensions to the perturbative 
theory to study the strong coupling structure of the bound states which will be the real 
challenge to this approach. 

Practically, even if we are successful in finding the light-front hadronic bound states 
in this direction, the theory is still not completely established. We have seen in Sec. Ill 
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that the matrix elements associated with various hadronic structmes arc defined on the 
canonical hght-front field operators in the light-front bound states. Since the bound states 
obtained in the above approach come from the use of the light-front effective QCD Hamil- 
tonian that contains many nontrivial counterterms for the severe light-front divergences 
that are allowed by the light-front power counting. Thus the matrix elements in terms 
of the canonical light-front field operators may not provide the correct physical informa- 
tion of hadrons due to the fact that many light-front physical operators associated with 
the obscrvablcs arc also dynamically dependent. Typical examples arc the transverse and 
longitudinal components of currents that have "bad" components ip- of fermions, and the 
operators in the hadronic structure function that belong to the high-twist contributions 
[94]. Therefore renormalization which introduces many noncanonical counterterms to the 
light-front Hamiltonian must also change the canonical structure of the operators used in 
the structure function definition. Obviously, the renormalization of physical operators be- 
sides the light-front Hamiltonian is another essential problem in the investigations of the 
light-front QCD bound states which has yet to be developed [156]. 



VIII. A BRIEF SUMMARY 



In this article, I have reviewed the basic formulation of the light-front field theory, 
the light-front expression of various strong interaction phenomena measured or extracted 
from various high energy experiments, and the possible understanding of these phenomena 
from QCD on the light-front. The essential problem for solving QCD on the light-front is to 
find the light-front hadronic bound states, since most of the interesting hadronic quantities, 
such as the parton distribution functions, parton fragmentation functions, hadronic form 
factors and coupling constants can naturally be formulated on the light-front in terms of 
simple matrix elements of light-front field operators in light-front hadronic bound states. 
Light-front field theory gives us a simple Schrodinger picture of bound state equations for 
relativistic particles, in which the bound states are simply a Fock states expansion on the 
trivial light-front vacuum. With such a picture, we may develop a reliable approach for 
addressing the relativistic bound states, as well as nonperturbative field theory for QCD. 
The main problem in this framework is reduced then, to constructing an effective QCD 
Hamiltonian on the light-front with a trivial vacuum so that it can reproduce the well- 
known phenomenology of hadrons on the low energy scale and can deduce perturbative 
QCD with asymptotic freedom at the high energy scale. A schematic approach to solving 
this problem has been discussed; further applications need to be explored. 
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APPENDIX A: TWO-COMPONENT FEYNMAN PERTURBATIVE 

THEORY 



For most field theorists Feynman's formulation is more convenient. In this appendix, 
I present the two-component light-front QCD Feynman theory. The two-component light- 
front QCD Lagrangian has a canonical form, 

£ = £o + Ant , (Al) 
where 



C-int — "Hint ■ 



(A2) 



The Feynman rules can be derived from the path integral formalism [158]. Here we use 
the procedure of 't Hooft and Veltman [159]. The result follows. The free quark and gluon 
propagators in two-component light-front QCD are 

^Oa(3{p) = ^2_^2 ' ^oib = ^iAh^ ■ (A3) 

The diagrammatic rules for various vertices are listed in Table III. The rest of the Feynman 
rules are the same as in the instant-form for boson and fermion theories. With these 
basic ingredients, the conventional Feynman approach to various physical processes can be 
directly discussed within the above two-component formulation of light-front QCD. 

It may be worth mentioning that one advantage of the two-component Feynman 
perturbation theory is that we can use the ML prescription and dimensional regulator to 
regulate the light-front infrared divergences and ultraviolet divergences in Feynman loop 
integrals, and may recover the multiplicative renormalizability of the theory, at least in the 
one-loop approximation, as shown in Ref. [60]. 
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TABLE III. Two-component Feynman Diagrammatic Rules 
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